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Milestone  results  have  been  achieved  in  the  design,  processing,  and  testing  of  dielectric  and  metallodielectric  composite  materials 
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Statement  of  the  problem  studied 

The  problem  studied  pertained  to  the  area  of  Highly-Controlled  Infrared  Dielectric  Emissivity.  In  response  to 
the  military  need  for  new  low-cost  IR  materials.  Novel  composite  materials,  such  as  photonic  crystals,  whose 
electromagnetic  properties  can  be  fully  engineered  and  manufactured  by  low-cost  techniques  were  needed  to  be 
developed. 


Summary  of  the  most  important  results 

The  MIT  Photonic  Crystal  HIDE  Materials  Engineering  Team  has  achieved  milestone  results  in  designing, 
processing  and  testing  dielectric  and  metallodielectric  composite  materials  for  control  of  IR  emissivity.  Our 
theoretical  analyses  have  discovered  that  an  omnidirectional  reflectivity  band  can  be  produced  with  high 
dielectric  contrast,  multilayer  structures,  as  discussed  in  detail  in  attachment  #1.  The  theory  was  applied  to 
design  and  fabricate  a  structure  comprised  of  alternating  Te  and  polymer  layers  that  were  produced  by  low  cost 
evaporation  and  spin  coating  processes,  respectively.  These  materials  exhibit  near  perfect,  omnidirectional 
reflectivity  in  the  8- 12  pm  range.  This  groundbreaking  result  has  received  international  acclaim  and  was 
published  in  Science  Magazine,  The  theoretical  and  experimental  details  of  this  1-D  omnidirectional  reflector 
are  presented  in  attachment  #  2.  We  have  also  developed  a  novel,  low-cost  process  for  the  periodic  placement  of 
micron-demensioned  spheres  and  clusters  of  metal  particles  in  two  dimensional  arrays,  as  detailed  in  attachment 
#3  and  #4.  This  process  is  the  first  example  of  stamping  and  filling  a  template  with  high  precision  at 
dimensions  relevant  to  metallodielectric  HIDE  materials.  A  sol-gel  technique  for  producing  multilayer  dielectric 
films  was  developed  and  led  to  the  fabrication  of  a  nearly  omnidirectional  reflector  corresponding  to 
1 ,55microns.  This  discovery  is  outlined  in  attachment  #5.  Our  team  used  the  omnidirectional  mirror  design 
criteria  to  develop  a  hollow  all-dielectric  waveguide  for  low-loss  transmission  of  electromagnetic  radiation.  The 
design  and  experimental  verification  of  this  novel  waveguide  Is  described  in  attachment  #6  This  in  turn  has  led 
to  the  development  of  a  novel  coaxial  optical  fiber  which  has  been  shown  theoretically  to  support  a  TEM  like 
mode  similar  to  the  one  existing  in  metallic  coaxial  cables.  This  discovery  was  published  in  Science  and  is 
described  in  detail  In  attachment  #7. 
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Omnidirectional  reflection  from  a  one-dimensional 

photonic  crystal 

Joshua  N.  Winn,  Yoel  Fink,  Shanhui  Fan,  and  J,  D,  Joannopoulos 

Massachusetts  Institute  of  Technology,  77  Massachusetts  Avenue,  Cambridge ,  Massachusetts  02139 
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We  demonstrate  that  one-dimensional  photonic  crystal  structures  (such  as  multilayer  films)  can  exhibit 
complete  reflection  of  radiation  in  a  given  frequency  range  for  all  incident  angles  and  polarizations.  We 
derive  a  general  criterion  for  this  behavior  that  does  not  require  materials  with  very  large  indices.  We 
perform  numerical  studies  that  illustrate  this  effect.  ©  1998  Optical  Society  of  America 
OCIS  codes:  230.4170,230.1480, 


Low-loss  periodic  dielectrics,  or  photonic  crystals,  allow 
the  propagation  of  light  to  be  controlled  in  otherwise 
difficult  or  impossible  ways.1'4  In  particular,  a  pho¬ 
tonic  crystal  can  be  a  perfect  mirror  for  light  from  any 
direction,  with  any  polarization,  within  a  specified  fre¬ 
quency  range.  It  is  natural  to  assume  that  a  necessary 
condition  for  such  omnidirectional  reflection  is  that  the 
crystal  exhibit  a  complete  three-dimensional  photonic 
bandgap,  that  is,  a  frequency  range  within  which  there 
are  no  propagating  solutions  of  Maxwell’s  equations. 
Here  we  report  that  this  assumption  is  false — in  fact  a 
one-dimensional  photonic  crystal  will  suffice.  We  in¬ 
troduce  a  general  criterion  for  omnidirectional  reflec¬ 
tion  for  all  polarizations  and  apply  it  to  the  case  of  a 
dielectric  multilayer  film.  Previous  attempts  to  attain 
high  reflectance  for  a  wide  range  of  incident  angles  in¬ 
volved  dielectric  films  with  high  indices  of  refraction, 
high  special  dispersion  properties,  or  multiple  contigu¬ 
ous  stacks  of  films.5'9 

A  one-dimensional  photonic  crystal  has  an  index 
of  refraction  that  is  periodic  in  the  y  coordinate  and 
consists  of  an  endlessly  repeating  stack  of  dielectric 
slabs,  which  alternate  in  thickness  from  dx  to  d%  and 
in  index  of  refraction  from  m  to  n%.  Incident  light 
can  be  either  s  polarized  (E  is  perpendicular  to  the 
plane  of  incidence)  or  p  polarized  (parallel).  Because 
the  medium  is  periodic  in  y  and  homogeneous  in  x  and 
z ,  the  electromagnetic  modes  can  be  characterized  by  a 
wave  vector  k,  with  ky  restricted  to  0  <  ky  <  tt fa.  We 
may  suppose  that  kz  =  (),&*>  0,  and  n2  >  nx  without 
loss  of  generality.  The  allowed  mode  frequencies  con 
for  each  choice  of  k  constitute  the  band  structure  of  the 
crystal.  The  continuous  functions  o>n(k),  for  each  n, 
are  the  photonic  bands. 

For  an  arbitrary  direction  of  propagation,  it  is  conve¬ 
nient  to  examine  the  projected  band  structure,  which  is 
shown  in  Fig.  1  for  a  quarter-wave  stack  with  nx  =  1 
and  t%2  ~  2.  To  make  this  plot  we  first  computed  the 
bands  con(kXr  ky)  for  the  structure,  using  a  numeri¬ 
cal  method  to  solve  Maxwell’s  equations  in  a  periodic 
medium.10  (In  fact,  for  the  special  case  of  a  multi¬ 
layer  film,  an  analytic  expression  for  the  dispersion 
relation  is  available.11)  Then,  for  each  value  of  kXf 
the  mode  frequencies  a>n  for  all  possible  values  of  ky 
were  plotted.  Thus  in  the  gray  regions  there  are  elec¬ 
tromagnetic  modes  for  some  value  of  kyj  whereas  in 
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the  white  regions  there  are  no  electromagnetic  modes, 
regardless  of  ky . 

One  obvious  feature  of  Fig.  1  is  that  there  is  no 
complete  bandgap.  For  any  frequency  there  exists 
some  electromagnetic  mode  with  that  frequency — the 
normal-incidence  bandgap  is  crossed  by  modes  with 
kz  >  0.  This  is  a  general  feature  of  one-dimensional 
photonic  crystals. 

However,  the  absence  of  a  complete  bandgap  does 
not  preclude  omnidirectional  reflection.  The  criterion 
is  not  that  there  be  no  propagating  states  within 
the  crystal;  rather,  the  criterion  is  that  there  be 
no  propagating  states  that  can  couple  to  an  incident 
propagating  wave.  As  we  argue  below,  the  latter 
criterion  is  equivalent  to  the  existence  of  a  frequency 
range  in  which  the  projected  band  structures  of  the 
crystal  and  the  ambient  medium  have  no  overlap. 

The  electromagnetic  modes  in  the  ambient  medium 
obey  a>  =  c(kx 2  +  ky2)mf  where  c  is  the  speed  of  light  in 
the  ambient  medium,  so  generally  w  >  ckx .  The  whole 
region  above  the  solid  diagonal  light  lines  at  ~  ckx  is 
filled  with  the  projected  bands  of  the  ambient  medium. 

If  a  semi-infinite  crystal  occupies  y  <  0  and  the 
ambient  medium  occupies  y  >  0,  the  system  is  no  longer 
periodic  in  the  y  direction  and  the  electromagnetic 
modes  of  the  system  can  no  longer  be  classified  by  a 
single  value  of  ky.  These  modes  must  be  written  as 


Parallel  wave  vector  kx  (2 nia) 

Fig.  1.  Projected  band  structure  for  a  quarter-wave  stack 
with  nx  —  1  and  n%  —  2.  Electromagnetic  modes  exist  only 
in  the  shaded  regions.  The  s -polarized  modes  are  plotted 
to  the  right  of  the  origin,  and  the  p-polarized  to  the  left. 
The  dark  lines  are  the  light  lines  to  —  ckx.  Frequencies 
are  reported  in  units  of  2irc/a , 
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a  weighted  sum  of  plane  waves  with  all  possible  ky. 
However,  kx  is  still  a  valid  symmetry  label.  The  angle 
of  incidence  0  upon  the  interface  at  y  =  0  is  related  to 
kx  by  co  sin  6  —  ckx . 

For  there  to  be  any  transmission  through  the  semi¬ 
infinite  crystal  at  a  particular  frequency,  there  must 
be  an  electromagnetic  mode  available  at  that  frequency 
that  is  extended  for  both  y  >  0  and  y  <  0.  Such  a 
mode  must  be  present  in  the  projected  photonic  band 
structures  of  both  the  crystal  and  the  ambient  medium. 
(The  only  states  that  could  be  present  in  the  semi¬ 
infinite  system  that  were  not  present  in  the  bulk 
system  are  surface  states,  which  decay  exponentially 
in  both  directions  away  from  the  surface  and  are 
therefore  irrelevant  to  the  transmission  of  an  external 
wave).  Therefore,  the  criterion  for  omnidirectional 
reflection  is  that  there  exist  a  frequency  zone  in  which 
the  projected  bands  of  the  crystal  have  no  states 
with  o)  >  ckx. 

In  Fig.  1,  the  lowest  two  p  bands  cross  at  a  point 
above  the  line  co  =  ckx ,  preventing  the  existence  of  such 
a  frequency  zone.  This  crossing  occurs  at  the  Brewster 
angle  9  b  ~  tan-1(«2/^i),  at  which  there  is  no  reflection 
of  p-polarized  waves  at  any  interface.  At  this  angle 
there  is  no  coupling  between  waves  with  ky  and  —ky>  a 
fact  that  permits  the  band  crossing  to  occur. 

This  difficulty  vanishes  when  we  lower  the  bands  of 
the  crystal  relative  to  those  of  the  ambient  medium  by 
raising  the  indices  of  refraction  of  the  dielectric  films. 
Figure  2  shows  the  projected  band  structure  for  the 
case  rii  =  1.7  and  n2  —  3.4.  In  this  case  there  is  a 
frequency  zone  in  which  the  projected  bands  of  the 
crystal  and  ambient  medium  do  not  overlap,  namely, 
from  the  filled  circle  (<oa/27rc  =  0.21)  to  the  open 
circle  {(oa/2iTC  =  0.27).  This  zone  is  bounded  above 
by  the  normal-incidence  bandgap  and  below  by  the 
intersection  of  the  top  of  the  first  gray  region  for 
p-polarized  waves  with  the  light  line. 

Between  the  frequencies  corresponding  to  the  filled 
and  open  circles  there  will  be  total  reflection  from 
any  incident  angle  for  either  polarization.  For  a  finite 
number  of  films  the  transmitted  light  will  diminish  ex¬ 
ponentially  with  the  number  of  films.  The  calculated 
transmission  spectra  for  a  finite  system  of  ten  films 
(five  periods)  are  plotted  in  Fig.  3  for  various  angles 
of  incidence.  The  calculations  were  performed  with 
transfer  matrices.12  The  stop  band  shifts  to  higher 
frequencies  with  more-oblique  angles,  but  there  is  a  re¬ 
gion  of  overlap  that  remains  Intact  for  all  angles. 

The  graphic  criterion  for  omnidirectional  reflection 
is  that  the  filled  circle  be  lower  than  the  open  circle 
(the  second  band  at  kx  —  0,  ky  =  7r /a).  Symbolically, 

&P\{kx  =  ~~^’ky  =  <  (op 2{kx  =  0,  ky  =  > 

(1) 

where  copn(kXt  ky)  Is  the  p-polarized  band  structure 
function  for  the  multilayer  film.  Note  that  the  left- 
hand  side  is  a  self-consistent  solution  for  frequency 
cOp\.  The  difference  between  these  two  frequencies  is 
the  range  of  omnidirectional  reflection. 


We  calculated  this  range  (when  It  exists)  for  a 
comprehensive  set  of  film  parameters.  Since  all  the 
mode  wavelengths  scale  linearly  with  dt  +  d2  =  a,  we 
need  consider  only  three  parameters  for  a  multilayer 
film:  Tii,  n2y  and  d\/a.  To  quantify  the  range  of  om¬ 
nidirectional  reflection  (aq,  &>2)  in  a  scale-independent 
manner,  we  report  the  range -midrange  ratio,  which  is 
defined  as  (&*2  -  )/V2(a>2  + 

For  each  choice  of  rii  and  n2/nly  there  is  a  value  of 
di/a  that  maximizes  the  range-midrange  ratio.  That 
choice  can  be  computed  numerically.  Figure  4  is  a 
contour  plot  of  the  ratio,  as  j%i  and  n2/n1  are  varied, 
for  the  maximizing  value  of  dtfa . 

An  approximate  analytic  expression  for  the  optimal 
zone  of  omnidirectional  reflection  can  be  derived: 


*  a  cos) 

A  co 

h 

fA-  2) 
A+2) 

|  a  cos| 

b 

IB-2\ 
B  +  2) 

2c  d\Tti  +  d2n2 

diyfni 2  - 

1  +  d2^Jri22 

'  ~  1 

(2) 


Fig.  2.  Projected  band  structure  for  a  quarter- wave  stack 
with  rii  —  1.7  and  n2  =  3.4,  with  the  same  conventions  as 
in  Fig.  1. 
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Fig.  3.  Calculated  transmission  spectra  for  a  quarter- 
wave  stack  of  ten  films  (ni  =  1.7,  n2  =  3.4)  for  three  angles 
of  incidence.  Solid  curves,  p-polarized  waves;  dashed 
curves,  s-polarized  waves.  The  overlapping  region  of  high 
reflectance  (>20  dB)  corresponds  to  the  region  between  the 
open  and  filled  circles  of  Fig.  2. 
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Smaller  index,  n1 

Fig.  4.  Range— midrange  ratio  for  omnidirectional  reflec¬ 
tion,  plotted  as  contours.  For  the  solid  contours  the  opti¬ 
mal  value  of  di/a  was  chosen.  The  dashed  curve  is  the 
0%  contour  for  the  case  of  a  quarter-wave  stack.  For  the 
general  case  of  an  ambient  medium  with  index  nQ  ^  1, 
the  abscissa  becomes 


where 

A=—  +  —,  B  =  n2^ni2  ~  1  |  niVn22  ~  1 

ni  n2  nxV«22  —  1  n-iyfn  i2  —  1 

(3) 

In  deriving  Eq.  (2)  we  assumed  that  the  optimal  film 
is  approximately  a  quarter-wave  stack.  Numerically 
we  find  this  to  be  an  excellent  approximation  for 
the  entire  range  of  parameters  depicted  in  Fig.  4;  the 
frequencies  as  predicted  by  this  approximation  are 
within  0.5%  of  the  exact  frequencies.  As  a  result  the 
optimization  of  d\j d  results  in  a  range -midrange  ratio 
very  close  to  that  which  results  from  a  quarter-wave 
stack  with  the  same  indices:  d\/a  =  n2/(n2  +  n\ ). 
In  Fig.  3,  the  0%  contour  for  quarter-wave  stacks 
is  plotted  as  a  dashed  curve,  which  is  very  close 
(always  within  2%  in  the  indices  of  refraction)  to  the 
numerically  optimized  contour. 

It  can  be  seen  from  Fig.  4  that,  for  omnidirectional 
reflection,  the  index  ratio  should  be  reasonably  high 
(>1.5)  and  the  indices  themselves  somewhat  higher 
(by  >1.5)  than  that  of  the  ambient  medium.  The 
former  condition  increases  the  band  splittings,  and 
the  latter  depresses  the  frequency  of  the  Brewster 
crossing.  An  increase  in  either  factor  can  partially 
compensate  for  the  other.  The  materials  should  also 
have  a  long  absorption  length  for  the  frequency  range 
of  interest,  especially  at  grazing  angles,  where  the  path 
length  of  the  reflected  light  along  the  crystal  surface 
is  long. 


Although  we  have  illustrated  our  arguments  by 
use  of  multilayer  films,  the  notions  in  this  Letter 
apply  generally  to  any  periodic  dielectric  function 
n{y).  What  is  required  is  the  existence  of  a  zone  of 
frequencies  in  which  the  projected  bands  of  the  crystal 
and  ambient  medium  have  no  overlap. 

However,  the  absence  of  a  complete  bandgap  does 
have  physical  consequences.  In  the  frequency  range 
of  omnidirectional  reflection  there  exist  propagating 
solutions  of  MaxwelTs  equations,  but  they  are  states 
with  a)  <  ckx  and  decrease  exponentially  away  from  the 
crystal  boundary.  If  such  a  state  were  launched  from 
within  the  crystal,  it  would  propagate  to  the  boundary 
and  reflect,  just  as  in  total  internal  reflection. 

Likewise,  although  it  might  be  arranged  that  the 
propagating  states  of  the  ambient  medium  do  not 
couple  to  the  propagating  states  of  the  crystal,  any 
evanescent  states  in  the  ambient  medium  will  couple  to 
them.  For  this  reason,  a  point  source  of  waves  placed 
very  close  (d  <  A)  to  the  crystal  surface  could  indeed 
couple  to  the  propagating  state  of  the  crystal.  Such 
restrictions,  however,  apply  only  to  a  point  source,  and 
one  can  easily  overcome  them  by  simply  adding  a  low- 
index  cladding  layer  to  separate  the  point  source  from 
the  film  surface. 
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A  Dielectric  Omnidirectional 
Reflector 

Yoel  Fink,  Joshua  N.  Winn,  Shanhui  Fan,  Chiping  Chen, 
Jurgen  Michel,  John  D.  Joannopouios,  Edwin  L.  Thomas* 


A  design  criterion  that  permits  truly  omnidirectional  reflectivity  for  all  polar¬ 
izations  of  incident  light  over  a  wide  selectable  range  of  frequencies  was  used 
in  fabricating  an  all-dielectric  omnidirectional  reflector  consisting  of  multilayer 
films.  The  reflector  was  simply  constructed  as  a  stack  of  nine  alternating 
micrometer-thick  layers  of  polystyrene  and  tellurium  and  demonstrates  om¬ 
nidirectional  reflection  over  the  wavelength  range  from  10  to  15  micrometers. 
Because  the  omnidirectionality  criterion  is  general,  it  can  be  used  to  design 
omnidirectional  reflectors  in  many  frequency  ranges  of  interest.  Potential  uses 
depend  on  the  geometry  of  the  system.  For  example,  coating  of  an  enclosure 
will  result  in  an  optical  cavity.  A  hollow  tube  will  produce  a  low-loss,  broadband 
waveguide,  whereas  a  planar  film  could  be  used  as  an  efficient  radiative  heat 
barrier  or  collector  in  thermoelectric  devices. 


Mirrors,  probably  the  most  prevalent  of 
optical  devices,  are  used  for  imaging  and 
solar  energy  collection  and  in  laser  cavities. 
One  can  distinguish  between  two  types  of 
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mirrors,  the  age-old  metallic  and  the  more 
recent  dielectric.  Metallic  mirrors  reflect 
light  over  a  broad  range  of  frequencies 
incident  from  arbitrary  angles  (that  is,  om¬ 
nidirectional  reflectance).  However,  at  in¬ 
frared  and  optical  frequencies,  a  few  per¬ 
cent  of  the  incident  power  is  typically  lost 
because  of  absorption.  Multilayer  dielectric 
mirrors  are  used  primarily  to  reflect  a  nar¬ 
row  range  of  frequencies  incident  from  a 
particular  angle  or  particular  angular  range. 
Unlike  their  metallic  counterparts,  dielec¬ 
tric  reflectors  can  be  extremely  low  loss. 
The  ability  to  reflect  light  of  arbitrary  angle 
of  incidence  for  all-dielectric  structures  has 
been  associated  with  the  existence  of  a 
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complete  photonic  band  gap  (i-J),  which 
can  exist  only  in  a  system  with  a  dielectric 
function  that  is  periodic  along  three  orthog¬ 
onal  directions.  In  fact,  a  recent  theoretical 
analysis  predicted  that  a  sufficient  condi¬ 
tion  for  the  achievement  of  omnidirectional 
reflection  in  a  periodic  system  with  an  in¬ 
terface  is  the  existence  of  an  overlapping 


band  gap  regime  in  phase  space  above  the 
light  cone  of  the  ambient  media  (4),  Now 
we  extend  the  theoretical  analysis  and  pro¬ 
vide  experimental  realization  of  a  multi¬ 
layer  omnidirectional  reflector  operable  in 
infrared  frequencies.  The  structure  is  made 
of  thin  layers  of  materials  with  different 
dielectric  constants  (polystyrene  and  tellu¬ 


rium)  and  combines  characteristic  features 
of  both  the  metallic  and  dielectric  mirrors. 
It  offers  metallic-like  omnidirectional  re¬ 
flectivity  together  with  frequency  selectiv¬ 
ity  and  low-loss  behavior  typical  of  multi¬ 
layer  dielectrics. 

We  consider  a  system  that  is  made  of  an 
array  of  alternating  dielectric  layers  cou¬ 
pled  to  a  homogeneous  medium,  character¬ 
ized  by  n0  (such  as  air  with  nQ  =  1),  at  the 
interface.  Electromagnetic  waves  are  inci¬ 
dent  upon  the  multilayer  film  from  the 
homogeneous  medium.  Although  such  a 
system  has  been  analyzed  extensively  in  the 
literature  (5-7),  the  possibility  of  omnidi¬ 
rectional  reflectivity  was  not  recognized 
until  recently.  The  generic  system  is  de¬ 
scribed  by  the  index  of  refraction  profile  in 
Fig.  1,  where  At  and  h2  are  the  layer  thick¬ 
ness  and  n%  and  n2  are  the  indices  of  refrac¬ 
tion  of  the  respective  layers.  The  incident 
wave  has  a  wave  vector  k  =  kxex  +  kyev  and 
a  frequency  of  u)  =  c\k\in0,  where  c  is  the 
speed  of  light  in  vacuum  and  ex  and  ey  are 
unit  vectors  in  the  x  and  y  directions,  re¬ 
spectively.  The  wave  vector  together  with 
the  normal  to  the  periodic  structure  defines 
a  mirror  plane  of  symmetry  that  allows  us 
to  distinguish  between  two  independent 
electromagnetic  modes:  transverse  electric 
(TE)  modes  and  transverse  magnetic  (TM) 
modes.  For  the  TE  mode,  the  electric  field 
is  perpendicular  to  the  plane,  as  is  the 
magnetic  field  for  the  TM  mode.  The  dis¬ 
tribution  of  the  electric  field  of  the  TE 
mode  (or  the  magnetic  field  in  the  TM 
mode)  in  a  particular  layer  within  the  strat¬ 
ified  structure  can  be  written  as  a  sum  of 
two  plane  waves  traveling  in  opposite  di¬ 
rections.  The  amplitudes  of  the  two  plane 
waves  in  a  particular  layer  a  of  one  cell  are 
related  to  the  amplitudes  in  the  same  layer 
of  an  adjacent  cell  by  a  unitary  2X2 
translation  matrix  (7). 

General  features  of  the  transport  prop¬ 
erties  of  the  finite  structure  can  be  under¬ 
stood  when  the  properties  of  the  infinite 
structure  are  elucidated.  In  a  structure  with 
an  infinite  number  of  layers,  translational 
symmetry  along  the  direction  perpendicular 
to  the  layers  leads  to  Bloch  wave  solutions 
of  the  form 

EK{x,y)  =  EK{x)eiK*e*#  (1) 

where  EK  (x,  v)  Is  a  field  component,  EK(x)  is 
periodic,  with  a  period  of  length  a,  and  K  is 
the  Bloch  wave  number  given  by 

K  =  ^ln(^  Tift/01) 

±  {j[Tr(t/“>)]2  -  lj  )  (2) 

Solutions  of  the  infinite  system  can  be 
propagating  or  evanescent,  corresponding 


Fig.  1.  Schematic  of  the  multi¬ 
layer  system  showing  the  layer 
parameters  {na  and  ha  are  the 
Index  of  refraction  and  thickness 
of  layer  a,  respectively),  the  in¬ 
cident  wave  vector  k,  and  the 
electromagnetic  mode  conven¬ 
tion.  E  and  B  are  the  electric  and 
magnetic  fields,  respectively. 
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Brewster  line 
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Parallel  wave  vector  ky  (2 n  /  a) 

Fig.  2.  (A)  Projected  band  structure  of  a  multilayer  film  with  the  light  line  and  Brewster  line, 
exhibiting  a  reflectivity  range  of  limited  angular  acceptance  with  n0  =  1,  /?1  =  2.2  and  nz  =  1.7 
and  a  thickness  ratio  of  h2/h1  =  2.2/ 1.7.  (B)  Projected  band  structure  of  a  multilayer  film  together 
with  the  light  line  and  Brewster  line,  showing  an  omnidirectional  reflectance  range  at  the  first  and 
second  harmonic.  Propagating  states,  light  gray;  evanescent  states,  white;  and  omnidirectional 
reflectance  range,  dark  gray.  The  film  parameters  are  /?1  =  4.6  and  n2  =  1.6  with  a  thickness  ratio 
of  hzih^  ~  1. 6/0.8.  These  parameters  are  similar  to  the  actual  polystyrene- tellurium  film 
parameters  measured  in  the  experiment. 
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to  real  or  imaginary  Bloch  wave  numbers, 
respectively.  The  solution  of  Eq.  2  defines 
the  band  structure  for  the  infinite  system, 
i ty(K,kv).  It  is  convenient  to  display  the  so¬ 
lutions  of  the  infinite  structure  by  project¬ 
ing  the  o>(Kyky)  function  onto  the  t o-ky 
plane;  Examples  of  such  projected  struc¬ 
tures  are  shown  in  Fig.  2,  A  and  B.  The 
light  gray  areas  highlight  phase  space 
where  K  is  strictly  real,  that  is,  regions  of 
propagating  states,  whereas  the  white  areas 
represent  regions  containing  evanescent 
states.  The  shape  of  the  projected  band 
structures  for  the  multilayer  film  can  be 
understood  intuitively.  At  ky  =  0,  the  band 
gap  for  waves  traveling  normal  to  the  lay¬ 
ers  is  recovered.  For  ky  >  0,  the  bands 
curve  upward  in  frequency.  As  k  — >  «  the 
modes  become  largely  confined  to  the  slabs 
with  the  high  index  of  refraction  and  do  not 
couple  between  layers  (and  are  therefore 
independent  of  £v). 

For  a  finite  structure,  the  translational 
symmetry  in  the  directions  parallel  to  the 
layers  is  preserved;  hence,  kv  remains  a 
conserved  quantity.  In  the  direction  perpen¬ 
dicular  to  the  layers,  the  translational  sym¬ 
metry  no  longer  exists.  Nevertheless,  the  K 
number,  as  defined  in  Eq.  2,  is  still  rele¬ 
vant,  because  it  is  determined  purely  by  the 
dielectric  and  structural  property  of  a  single 
bilayer.  In  regions  where  K  is  imaginary, 
the  electromagnetic  field  is  strongly  atten¬ 
uated.  As  the  number  of  layers  is  increased, 
the  transmission  coefficient  decreases  ex¬ 
ponentially,  whereas  the  reflectivity  ap¬ 
proaches  unity. 

Because  we  are  primarily  interested  in 
waves  originating  from  the  homogeneous 
medium  external  to  the  periodic  structure, 
we  will  focus  only  on  the  portion  of  phase 
space  lying  above  the  light  line.  Waves 
originating  from  the  homogeneous  medium 
satisfy  the  condition  o>  >:  ckjn0y  where  nQ  is 
the  refractive  index  of  the  homogeneous 
medium,  and  therefore  they  must  reside 
above  the  light  line.  States  of  the  homoge¬ 
neous  medium  with  k  —  0  are  normal 
incident,  and  those  lying  on  the  =  ck Jn0 
line  with  kx  =  0  are  incident  at  an  angle  of 
90°. 

The  states  in  Fig,  2A  that  are  lying  in 
the  restricted  phase  space  defined  by  the 
light  line  and  that  have  a  (w,  ky)  corre¬ 
sponding  to  the  propagating  solutions  (gray 
areas)  of  the  crystal  can  propagate  in  both 
the  homogeneous  medium  and  the  struc¬ 
ture.  These  waves  will  partially  or  en¬ 
tirely  transmit  through  the  film.  Those 
states  with  (<*>,  &v)  in  the  evanescent  regions 
(white  areas)  can  propagate  in  the  homoge¬ 
neous  medium  but  will  decay  in  the  crys¬ 
tal — waves  corresponding  to  this  portion 
of  phase  space  will  be  reflected  off  the 
slructure. 
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The  multilayer  system  leading  to  Fig. 
2A  represents  a  structure  with  a  limited 
reflectivity  cone  because  for  any  frequency 
one  can  always  find  a  kv  vector  for  which  a 
wave  at  that  frequency  can  propagate  in  the 
crystal  and  hence  transmit  through  the  film. 
For  example,  a  wave  with  w  =  0.285  X 
I'Kda  (dashed  horizontal  line  in  Fig.  2A) 
will  be  reflected  for  a  range  of  ky  values 
ranging  from  0  (normal  incidence)  to 
0.285  X  2Wcr  (90°  incidence)  in  the  TE 
mode,  whereas  in  the  TM  mode  it  begins  to 
transmit  at  a  value  of  ky  =  0.187  X  2n7u 
(—41°  incidence).  The  necessary  and  suffi¬ 
cient  criterion  (8)  for  omnidirectional  re¬ 
flectivity  at  a  given  frequency  is  that  no 
transmitting  state  of  the  structure  exists 
inside  the  light  cone;  this  criterion  is  satis¬ 
fied  by  frequency  ranges  marked  in  dark 
gray  in  Fig.  2B.  In  fact,  the  system  leading 
to  Fig.  2B  exhibits  two  omnidirectional 
reflectivity  ranges. 

The  omnidirectional  range  is  defined 
from  above  by  the  normal  incidence  band 
edge  &h(kx  —  tt fa,  ky  =  0)  (point  a  in  Fig. 
2B)  and  from  below  by  the  intersection  of 
the  top  of  the  TM  allowed  band  edge  with 
the  light  line  iaf(kx  =  rr/a,  ky  =  ru/c)  (point 
b  in  Fig.  2B). 

The  exact  expression  for  the  band  edges  is 
^Y~cos  (kjl)hl  +  kxmh2) 

+  — 2 —  co s(kxinh}  -  kj2)h2)  +  1=0,  (3) 

where  kxi<x>  =  V(oj njcf  ~  J?  (a  =  1,2) 
and 


Ukxm 
2\^  + 
»,2k,'2) 
n2%{"  + 


kjn 
k™ 

n2k}X) 


TE 


TM 


(4) 


A  dimensionless  parameter  used  to  quantify  the 
extent  of  the  omnidirectional  range  is  the  range 
to  midrange  ratio  defined  as  {b)h  -  + 

w,).  Figure  3  is  a  plot  of  this  ratio  as  a  taction  of 
njn}  and  nxfnQ,  where  and  u>/  are  determined 
by  solutions  of  Eq.  3  with  quarter  wave  layer 
thickness.  The  contours  in  this  figure  represent 
various  equiomnidireetional  ranges  for  different 
material  index  parameters  and  could  be  useful  for 
design  purposes. 

It  may  also  be  useful  to  have  an  approxi¬ 
mate  analytical  expression  for  the  extent  of 
the  gap.  This  can  be  obtained  by  setting 
cos (V%  -  kj2>h2)  s  1  in  Eq.  3.  We  find 
that  for  a  given  incident  angle  0O,  the  approx¬ 
imate  width  in  frequency  is 

Au>(0o) 


2c 

hi  ^n,2  -  iio2sm20o  +  h2  yjn2  ~  n02s\n2B0 


(5) 


At  normal  incidence,  there  is  no  distinction 
between  TM  and  TE  modes.  At  increasingly 
oblique  angles,  the  gap  of  the  TE  mode  in¬ 
creases,  whereas  the  gap  of  the  TM  mode 
decreases.  In  addition,  the  center  of  the  gap 
shifts  to  higher  frequencies.  Therefore,  the 


100 


frequency  range  of  omnidirectional  reflection, 
plotted  as  contours.  Here,  the  layers  were  set 
to  quarter  wave  thickness  and  >  n2.  The 
ratio  for  our  materials  is  about  45%  (n%/n2  — 

2.875  and  nzInQ  =  1.6}.  It  is  located  at  the 
intersection  of  the  dashed  lines  (black  dot).  Fig.  4  (right).  Calculated  (solid  line)  and  measured 
(dashed  line)  reflectance  (in  percent)  as  a  function  of  wavelength  for  TN  and  TE  modes  at  normal, 
45°,  and  80°  angles  of  incidence,  showing  an  omnidirectional  reflectivity  band. 
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criterion  for  the  existence  of  omnidirectional 
reflectivity  can  be  restated  as  the  occurrence 
of  a  frequency  overlap  between  the  gap  at 
normal  incidence  and  the  gap  of  the  TM 
mode  at  90°.  Analytical  expressions  for  the 
range  to  midrange  ratio  can  be  obtained  by 
setting 


2c 


to*  = 


h2n2  + 
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2c 
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h2yfn22  -  nQ2  4-  hx  ylnf"-  n02 
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„  2  CT" 

^2  Y^l 


«i2V”22  -  «02  +  «22\/*l2  -  «02 


(6b) 

Moreover,  the  maximum  range  width  is  at¬ 
tained  for  thickness  values  that  are  not  equal 
to  the  quarter  wave  stack  although  the  in¬ 
crease  in  band  width  gained  by  deviating 
from  the  quarter  wave  stack  is  typically  only 
a  few  percent  (4). 

In  general,  the  TM  mode  defines  the  low¬ 
er  frequency  edge  of  the  omnidirectional 
range.  An  example  can  be  seen  in  Fig.  2B  for 
a  particular  choice  of  the  indices  of  refrac¬ 
tion.  This  can  be  proven  by  showing  that 

dix)  I  £)co  I 

dky  I  tm  dkv  I  te  ^  ^ 

in  the  region  that  resides  inside  the  light 
line.  The  physical  reason  for  Eq.  7  lies  in 
the  vectorial  nature  of  the  electric  field.  In 
the  upper  portion  of  the  first  band,  the 
electric  field  concentrates  its  energy  in  the 
high  dielectric  regions.  Away  from  normal 
incidence,  the  electric  field  in  the  TM  mode 
has  a  component  in  the  direction  of  period¬ 
icity,  and  this  component  forces  a  larger 
portion  of  the  electric  field  into  the  low 
dielectric  regions.  The  group  velocity  of 
the  TM  mode  is  therefore  enhanced.  In 
contrast,  the  electric  field  of  the  TE  mode 
is  always  perpendicular  to  the  direction  of 
periodicity  and  can  concentrate  its  energy 
primarily  in  the  high  dielectric  region. 

A  polystyrene-tellurium  (PS-Te)  materi¬ 
als  system  was  chosen  to  demonstrate  omni¬ 
directional  reflectivity.  Tellurium  has  a  high 
index  of  refraction  and  low  loss  characteris¬ 
tics  in  the  frequency  range  of  interest.  In 
addition,  its  relatively  low  latent  heat  of  con¬ 
densation  together  with  the  high  glass  transi- 


Tabte  1.  Penetration  depth  (£)  at  different  angles 
of  incidence  for  the  TE  and  TM  modes. 


Angle  of  incidence 
(degrees) 

(M-m) 

Ste  (m-™) 

0 

2,51 

2.51 

45 

3.05 

2.43 

80 

4.60 

2.39 
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tion  temperature  of  the  PS  minimizes  diffu¬ 
sion  of  Te  into  the  polymer  layer.  The  choice 
of  PS,  which  has  a  series  of  absorption  peaks 
in  the  measurement  range  (V),  demonstrates 
the  competition  between  reflectivity  and  ab¬ 
sorption  that  occurs  when  an  absorption  peak 
is  located  in  the  evanescent  state  region.  The 
Te  (0.8  jxm)  and  PS  (1.65  pm)  films  were 
deposited  (10)  sequentially  to  create  a  nine- 
layer  film  (11), 

The  optical  response  of  this  particular  mul¬ 
tilayer  film  was  designed  to  have  a  high  reflec¬ 
tivity  region  in  the  10-  to  15-pm  range  for  any 
angle  of  incidence  (in  the  experiment,  we  mea¬ 
sure  from  0°  to  80°).  The  optical  response  at 
oblique  angles  of  incidence  was  measured  with 
a  Fourier  Transform  Infrared  Spectrometer 
(Nicolet  860)  fitted  with  a  polarizer  (ZnS;  Spec- 
IraTech)  and  an  angular  reflectivity  stage 
(VeeMax;  SpectraTech).  At  normal  incidence, 
the  reflectivity  was  measured  with  a  Nicolet 
Infrared  Microscope.  A  freshly  evaporated  alu¬ 
minum  mirror  was  used  as  a  background  for  the 
reflectance  measurements. 

Good  agreement  between  the  calculated 
(12)  and  measured  reflectance  spectra  at 
normal,  45°,  and  80°  incidence  for  the  TM 
and  TE  modes  is  shown  in  Fig.  4.  The 
regimes  of  high  reflectivity  at  the  different 
angles  of  incidence  overlap,  thus  forming  a 
reflective  range  of  frequencies  for  light  of 
any  angle  of  incidence.  The  frequency  lo¬ 
cation  of  the  omnidirectional  range  is  de¬ 
termined  by  the  layer  thickness  and  can  be 
tuned  to  meet  specifications.  The  range  is 
calculated  from  Eq.  6  to  be  5.6  pm,  and  the 
center  wavelength  is  12.4  pm,  correspond¬ 
ing  to  a  45%  range  to  midrange  ratio  shown 
in  dashed  lines  in  Fig.  3  for  the  experimen¬ 
tal  index  of  refraction  parameters.  These 
values  are  in  agreement  with  the  measured 
data.  The  calculations  are  for  lossless  me¬ 
dia  and  therefore  do  not  predict  the  PS 
absorption  band  at  —13  and  14  pm.  The  PS 
absorption  peak  is  seen  to  increase  at  larger 
angles  of  incidence  for  the  TM  mode  and  to 
decrease  for  the  TE  mode.  The  physical 
basis  for  these  phenomena  lies  in  the  rela¬ 
tion  between  the  penetration  depth  and  the 
amount  of  absorption.  The  penetration  depth  is 
|  «  Im(l/AT),  where  K  is  the  Bloch  wave  num¬ 
ber.  It  can  be  shown  that  £  is  a  monotonlcally 
increasing  function  of  the  incident  angle  for  the 
TM  mode  of  an  omnidirectional  reflector  and  is 
relatively  constant  for  the  TE  mode.  Thus,  the 
TM  mode  penetrates  deeper  into  the  structure  at 
increasing  angles  of  incidence  (Table  1)  and  is 
more  readily  absorbed.  The  magnitude  of  the 
imaginary  part  of  the  Bloch  wave  number  for  a 
mode  lying  in  the  gap  is  related  to  its  distance 
from  the  band  edges.  This  distance  increases 
in  the  TE  mode  because  of  the  widening  of 
the  gap  at  Increasing  angles  of  incidence  and 
decreases  in  the  TM  mode  because  of  the 
shrinking  of  the  gap. 


The  PS-Te  structure  does  not  have  a  com¬ 
plete  photonic  band  gap.  Its  omnidirectional 
reflectivity  is  due  instead  to  the  restricted 
phase  space  available  to  the  propagating 
states  of  the  system.  The  materials  and  pro¬ 
cesses  were  chosen  for  their  low  cost  and 
applicability  to  large  area  coverage.  The  pos¬ 
sibility  of  achieving  omnidirectional  reflec¬ 
tivity  itself  is  not  associated  with  any  partic¬ 
ular  choice  of  materials  and  can  be  applied  to 
many  wavelengths  of  interest.  Our  structure 
offers  metallic-like  omnidirectional  reflectiv¬ 
ity  for  a  wide  range  of  frequencies  and  at  the 
same  time  is  of  low  loss.  In  addition,  it  allows 
the  flexibility  of  frequency  selection. 
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Summary 

Researchers  at  the  Ceramics  Processing  Research  Lab  are  investigating  non-conventional 
methods  of  particle  arrangement  in  an  effort  to  create  a  photonic  bandgap  material.  Experiments 
using  silica  spheres  and  copper  particles  placed  in  etched  silicon  substrates  have  been  conducted 
which  demonstrate  the  feasibility  of  precisely  arranging  particles.  A  micromolding  technique 
developed  at  the  CPRL  shows  promise  in  being  able  to  precisely  embed  particles  in  a  variety  of 
materials. 
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Figure  la:  Etched  silicon  substrate  used  as  the  master  Figure  lb:  Glass  spheres  precisely  ordered  into  hole 
tool  to  form  ordered  arrays  of  holes  in  cheaper  substrates  array.  The  glass  spheres  are  readily  obtainable  in  very 
such  as  plastic  or  wax.  narrow  size  distributions  and  are  used  in  this  case  to 

demonstrate  the  feasibility  of  the  process.  Metallized 
glass  spheres  or  narrow  size  metal  spheres  will  be  used 
in  the  final  composite 
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Experiments  with  monosized  glass  spheres 

Experiments  filling  the  holes  with  uniformly  sized  glass  spheres  have  been  very 
successful.  A  solution  of  ethanol  and  1.58  |xm  silica  standard  spheres  is  applied  to  a  silicon 
substrate  with  holes  approximately  2  pm  in  diameter  3.5  pm  apart.  The  substrate  is  sitting  on  a 
plate  in  a  sonic  bath  to  shake  the  particles  into  the  holes.  The  particles  settle  into  the  holes,  and 
the  excess  is  wiped  clean  after  the  ethanol  evaporates.  Several  applications  result  in  complete 
filing  of  the  matrix.  Figure  la  is  an  SEM  photograph  of  a  silicon  substrate,  which  can  be  used  as 
a  tool  to  mass  produce  matrices  out  of  other  materials  such  as  wax  or  plastic.  Figure  lb  shows  a 
completely  filled  array  of  glass  spheres,  demonstrating  the  feasibility  of  the  concept  of  precisely 
arranging  particles  in  an  array.  These  encouraging  results  show  promise  for  development  of  a 
system  of  reusable  tooling  able  to  capture  and  precisely  arrange  particles. 

Experiments  with  wide  size  distribution  metal  particles 

Experiments  in  filling  arrays  were  also  conducted  using  copper  powder.  The  powder 
used  has  a  wide  size  distribution,  so  many  holes  contain  several  particles  while  others  only 
contain  one  (Figure  2).  Improvements  in  powder  processing  are  required  to  yield  a  more  uniform 
ideal  particle  size  of  two  microns  in  order  to  perfectly  fill  the  holes.  Ideally,  a  uniform  particle 
size  material  would  fill  the  lattice  perfectly,  as  was  demonstrated  with  the  silica  spheres. 
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Figure  2:  Silicon  Substrate  Filled  With  1  pm  Classified  Copper  Particles 


Polymer  substrates 

A  micromolding  process  can  be  used  to  transfer  the  pattern  of  holes  to  other  materials. 
Preliminary  experiments  using  paraffin  have  been  conducted,  although  results  indicate  that 
paraffin  is  too  soft  a  media  to  apply  the  particles  too.  Wiping  of  the  substrate  between 
applications  damages  the  mold  (Figure  3).  Other  materials  such  as  epoxies  and  polyimides  are 
being  investigated  as  suitable  molding  materials.  The  ultimate  goal  is  to  create  a  thin  film 
containing  one  layer  of  the  lattice,  and  then  to  stack  several  films  together  to  form  a  diamond 
lattice  array. 
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Future  work 

Future  work  will  consist  of  investigating  different  materials  in  which  to  place  metallic 
particles,  searching  for  a  reliable  source  or  method  of  production  of  uniform  particle  size 
particles,  and  examining  the  challenges  involved  in  creating  several  layers  of  arranged  particles 
to  form  a  diamond  lattice.  Investigation  into  the  feasibility  of  metallizing  readily  available 
uniformly  sized  particles,  such  as  the  glass  spheres  used  in  this  work,  will  also  be  conducted. 
Current  processing  technology  is  unable  to  provide  monosized  metallic  particles,  and  is 
considered  a  barrier  to  current  research. 
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Metallodielectric  Photonic  Crystals 

Collaborators:  Shanhui  Fan,  Pierre  Villenueve,  Joshua  Weitz,  Matt  Rosenthal, 
Malinda  Tupper,  Michael  Cima  and  John  Joannopoulos 

Using  a  finite  difference  time-domain  method  we  studied  the  band  structure  of 
three-dimensional  metallodielectric  photonic  crystals.  Complete  gaps  are  found  in  a 
diamond  lattice  of  isolated  metal  spheres  as  well  as  for  metal  spheres  embedded  in  a 

dielectric  material  like  Teflon  (e  =  2.1).  The  gap  can  be  as  large  as  60%  when  the  radius  of 
spheres  is  approximately  21%  of  the  cubic  unit  cell  size. 

Such  optimal  gaps  require  that  spheres  overlap  across  layers.  We  have,  however, 
found  a  configuration  which  retains  a  large  omnidirectional  gap  and  can  be  fabricated  in  a 
layer-by-layer  fashion.  The  band  structure  of  this  configuration  using  metal  spheres  with 
radii  17.7%  of  the  cubic  unit  cell  size  is  presented  in  the  enclosed  figure.  The  fabrication 
design  is  currently  being  implemented  by  Matt  Rosenthal,  Malinda  Tupper  and  Michael 
Cima. 


The  unit  cell  size  can  be  scaled  depending  on  the  wavelength  of  interest  For 
example,  a  diamond  lattice  with  4.95  |im  unit  cell  size  and  non-overlapping  metal  spheres 

will  have  a  20%  gap  with  a  midgap  frequency  of  8.8  pm.  To  accommodate  layer  by  layer 

growth  the  particles  should  be  mono-dispersed  with  diameter  1.75  |im.  Transmission 
calculations  predict  a  6  dB  rejection  per  layer  for  isolated  metal  spheres. 

We  are  also  developing  a  finite  difference  time  domain  code  that  will  model  the 
dispersion  relationships  of  more  complicated  materials.  Using  this  method  we  have 
successfully  calculated  the  dispersion  curves  of  polariton  and  drude  metals.  Further  testing 
is  necessary  to  confirm  that  the  simulations  correctly  model  the  time  decay  of  propagating 
fields. 


Frequency  (c/a) 
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Band  Structure  for  a  Metallodielectric  Photonic  Crystal 
Arranged  in  a  Diamond  Lattice  with  Non-Overlapping  Metal  Spheres 
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Thin  films  of  Si02  and  Ti02  were  used  to  fabricate  one-dimensional  photonic  crystal  devices  using 
the  sol-gel  method:  an  omnidirectional  reflector  and  microcavity  resonator.  The  reflector  consisted 
of  six  Si02/Ti02  bilayers,  designed  with  a  stopband  in  the  near  infrared.  Reflectivity  over  an 
incident  angle  range  of  0°— 80°  showed  an  omnidirectional  band  of  70  nm,  which  agrees  with 
theoretical  predictions  for  this  materials  system.  The  microcavity  resonator  consisted  of  a  Ti02 
Fabry- Perot  cavity  sandwiched  between  two  Si02/Ti02  mirrors  of  three  bilayers  each.  We  have 
fabricated  a  microcavity  with  resonance  at  X  =  1 500  nm  and  achieved  a  quality  factor  of  Q  =  35.  We 
measured  a  resonance  frequency  modulation  with  a  change  in  incident  angle  of  light  and  defect  layer 
thickness.  ©  1999  American  Institute  of  Physics.  [§ 0003-6951(99)03850-4] 


The  fabrication  techniques  typically  used  for  solid  state 
microphotonic  components  are  well  suited  for  deposition 
onto  wafer-sized  surface  areas.  When  deposition  onto  larger 
areas  is  required,  other  techniques  that  are  better  suited  for 
macrophotonic  applications  may  provide  practical  advan¬ 
tages.  An  example  of  such  a  technique  is  the  sol-gel 
method.1  The  sol-gel  method  is  an  inexpensive  wet- 
chemistry  process  that  has  been  used  in  micro-  and  macro¬ 
photonic  realms  to  deposit  a  variety  of  films  onto  various 
substrates  as  antireflection  (AR)  coatings,1  selective  filters,2 
Er-doped  optical  waveguides,3  and  Bragg  grating 
waveguides.4  The  sol-gel  method  offers  compositional  flex¬ 
ibility  over  a  range  of  oxides,  piezoelectrics,  and  ferroelee- 
trics. 

This  variety  of  materials  makes  it  possible  to  explore  the 
use  of  composite  sol-gel  structures  as  photonic  band-gap  de¬ 
vices.  As  optical  analogs  to  semiconductors,  photonic  crys¬ 
tals  possess  a  periodic  refractive  index  that  can  give  rise  to 
bands  of  forbidden  photon  propagation  states  in  the 
material.5  With  appropriate  design,  the  presence  of  a  defect 
in  the  periodic  structure  can  introduce  localized  photonic 
states  within  the  photonic  band  gap.  Examples  of  one¬ 
dimensional  photonic  crystals  are  an  omnidirectional 
reflector6  and  a  waveguide  microcavity.7  Interference  filters 
and  microcavities  fabricated  from  solid-state  materials  have 
been  well-studied  and  characterized. 

In  this  work,  we  demonstrate  the  feasibility  of  using  the 
sol-gel  method  to  fabricate  one-dimensional  photonic  crys¬ 
tals  with  and  without  a  defect — a  microcavity  and  omnidi¬ 
rectional  reflector,  respectively — for  macrophotonic  applica¬ 
tions.  We  choose  Si02  and  Ti02(nsi02=  1.47,  ^xi(>2=2,36, 
measured  at  632,8  nm)  as  a  materials  system  with  index 
contrast  sufficiently  large  for  omnidirectional  reflectance.6 
The  devices  that  we  fabricate,  because  they  are  composed  of 
oxides,  may  be  well-suited  for  application  in  environments 
that  see  elevated  temperatures,  corrosion,  and  radiation. 


^Electronic  mail:  lckim@mit.edu 


In  general,  the  sol-gel  method  begins  with  a  step- 
polymerization  process  wherein  an  organometallie  precursor 
is  hydrolyzed  to  form  a  metal-oxide  network  gel.  After  depo¬ 
sition,  typically  via  dip  or  spin  coating,  the  solution  is  dried 
and  fired  to  result  in  a  final  film  with  porosity  as  low  as  a  few 
percent  depending  on  the  heat  treatment. 

Solutions  of  Si02  and  Ti02  were  synthesized  using  reci¬ 
pes  requiring  minimal  chemistry  apparatus,  following  the 
methods  given  in  Refs,  8  and  9,  respectively.  All  chemical 
preparation  was  done  in  a  glovebox  filled  with  Ar  atmo¬ 
sphere.  Each  solution  was  magnetically  stirred  overnight  be¬ 
fore  use. 

To  deposit  a  thin  film,  several  milliliters  of  solution  were 
spin  cast  onto  a  cleaned  Si  wafer  piece.  The  sample  was  then 
placed  in  a  conventional  tube  furnace,  dried  at  100°C  for  2 
min,  and  fired  at  700  °C  for  15  min.  The  film  was  allowed  to 
cool  to  200  °C  for  10  min  before  removal  from  the  furnace. 
Films  of  Si02  and  Ti02  received  identical  heat  treatments. 

For  calibration,  film  thickness  and  refractive  index  were 
measured  using  a  Gaertner  Ellipsometer  operating  at  632.8 
nm.  Actual  thicknesses  were  obtained  via  cross-sectional 
transmission  electron  microscopy  (XTEM).  Reflectance 
measurements  were  done  using  a  Fourier  transform  infrared 
spectrometer  (Nicolet  860;  -2  mm2  beam  spot  size  in  nor¬ 
mal  incidence)  with  angle-dependence  stage  (SpectraTech 
Veemax;  ^20  mm2  spot  size  for  30°-80°  off-normal)  and 
ZnS  polarizer  (SpectraTech).  Each  measurement  was  refer¬ 
enced  against  an  Au  background. 

The  interference  filter  is  composed  of  six  bi layers  of 
Si02  and  Ti02  quarter-wave  films.  The  microcavity  is  com¬ 
posed  of  a  half-wave  Ti02  “defect”  layer  sandwiched  be¬ 
tween  mirrors  consisting  of  three  bilayers  each.  The  XTEM 
micrographs  of  these  devices  show  good  thickness  unifor¬ 
mity  for  individual  films  as  well  as  the  presence  of  interface 
and  film  defects  (Fig.  1). 

To  reduce  the  thermal  budget  on  the  interference  filter 
and  microcavity  devices,  a  minimum  temperature  for  densi- 
fieation  was  desired.  Single  films  of  Si02  and  TiG2  fired  at 
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FIG.  2.  Angle-dependent  reflectance  spectra  of  the  omnidirectional  reflector 
for  TE  and  TM  modes  (dotted  lines  show  calculated  spectra).  Omnidirec¬ 
tional  reflection  bands  are  marked  by  the  dotted  vertical  lines. 


w  mmrnmmmmmm 

500  nm 

FIG.  1 .  Cross-sectional  TEM  micrographs  of  devices  fabricated  by  the  sol- 
gel  method:  (a)  omnidirectional  reflector  and  (b)  microcavity  resonator.  The 
micrographs  detail  films  in  the  interior  of  the  devices. 

500  °C  contained  macroscopic  defects  thought  to  be  pores 
and  cracks  as  shown  by  optical  microscopy,  which  suggested 
incomplete  densification.  At  700  °C,  these  defects  were  not 
observed.  Moreover,  interference  filters  fabricated  using 
500  °C  anneals  showed  considerably  poorer  reflectance  spec¬ 
tra  than  those  fabricated  at  700  °C,  presumably  due  to  scat¬ 
tering  from  the  defects  and  thickness  irregularities.  There¬ 
fore,  700  °C  densification  anneals  were  used  for  all  devices. 

The  XTEM  micrograph  illustrates  that  the  spin-coated 
films  have  good  thickness  uniformity  and  that  the  deposition 
process  is  repeatable.  There  is  no  stratification  of  film  den¬ 
sity  due  to  the  accumulation  of  heat  treatments  for  initial 
films  deposited  in  the  multilayer,  which  would  lead  to  an 
increase  in  refractive  index8  and  appear  as  a  nonuniformity 
in  film  thickness  through  a  cross  section.  It  can  be  concluded 
that  the  15  min  firing  at  700  °C  is  sufficient  for  full  densifi¬ 
cation  of  these  films. 

The  angle-dependent  reflectance  plots  for  the  interfer¬ 
ence  filter  are  shown  in  Fig.  2.  For  the  transverse  electric 
(TE)  and  transverse  magnetic  (TM)  polarizations,  the  reflec¬ 
tance  bands  shift  to  shorter  wavelengths  with  increasing  in¬ 
cident  angle.  The  bandwidths,  however,  behave  differently: 
the  TE  bandwidths  increase  while  the  TM  bandwidths  de¬ 
crease  with  increasing  incident  angle.  The  reflectance  behav¬ 
ior  agrees  with  theoretical  dispersion  curves.6,10  Using  mea¬ 
sured  thicknesses  and  refractive  indices  for  Si02  and  Ti02, 
we  calculated  theoretical  angle-dependent  reflectance  spectra 
with  the  transfer  matrix  method.10  In  general,  there  is  good 
agreement  between  the  calculated  and  measured  spectra. 

Figure  2  shows  that  there  is  an  overlapping  reflectance 
band  common  to  all  spectra  taken  between  the  incident 
angles  of  0°— 80°;  this  is  the  omnidirectional  reflectance 
band.  For  TE  polarization,  there  is  a  300  nm  omnidirectional 
reflection  band.  For  TM  polarization  there  is  a  70  nm  band. 
The  superposition  of  TE  and  TM  polarizations  results  in  a  70 

nm  band  of  omnidirectional  reflectance  for  unpolarized  inci- 
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dent  light.  These  results  are  in  good  agreement  with  the  pre¬ 
dicted  omnidirectional  band6  for  this  materials  system,  which 
is  approximately  70  nm. 

The  ideal  design  for  highest  reflectance  and  sharpest  re¬ 
flectance  band  edges  in  an  interference  filter  is  a  quarter- 
wave  stack  because  Bragg  interference  among  propagating 
waves  within  the  structure  is  maximized.  When  film  thick¬ 
nesses  are  detuned  from  the  quarter-wave  optical  thick¬ 
nesses,  as  is  the  case  here,  Bragg  interference  is  weakened 
and  the  filter  exhibits  a  high  reflectance  range  (HRR)  sensi¬ 
tive  to  the  average  optical  thickness  of  the  layers  in  the 
filter.11  Detuning  can  relax  thickness  constraints  on  indi¬ 
vidual  layers  while  maintaining  high  reflectance  over  a  wider 
range  than  achievable  using  quarter-wave  thickness.  The 
tradeoff  is  that  there  is  a  decrease  in  the  maximum  reflec¬ 
tance  as  well  as  a  reduction  in  sharpness  of  the  reflectance 
band  edges.  Applications  requiring  a  wide  HRR  and  process 
flexibility  stand  to  benefit  from  detuning. 

The  reflectance  spectra  of  a  microcavity,  on  the  other 
hand,  is  less  forgiving  to  thickness  and  index  variations.  In  a 
microcavity,  the  defect  layer  accommodates  a  resonant  opti¬ 
cal  mode  whose  linewidth  is  determined  by  the  reflectances 
of  the  adjacent  mirrors.  Since  the  cavity  quality  factor  Q  is 
sensitive  to  the  mirror  reflectances  Rx  and  R2  as  Q 
|  In  |  2 1  1  (f°r  a  lossless  Fabry-Perot  cavity),  detuning 
in  the  mirror  layers  decreases  Q  by  increasing  leakage  of  the 
confined  mode.  Thus  by  measuring  Q,  we  obtain  a  direct 
indication  of  the  thickness  and  index  accuracy  of  films  de¬ 
posited  by  the  sol-gel  method. 

Angle-dependent  reflectance  plots  for  the  microcavity 
are  shown  in  Fig.  3.  In  normal  incidence,  the  cavity  reso¬ 
nance  occurs  at  X=1500  nm  and  a  highest  quality  factor  of 
2  =  35  was  measured.  By  comparison,  solid-state  vacuum 
deposition  techniques  can  be  used  to  fabricate  microcavities 
with  Q—  10000.  Figure  1(b)  illustrates  that  film  thicknesses 
for  the  microcavity  showed  more  variations  than  seen  in  the 
interference  filter.  The  low  Q  indicates  that  the  deposition 
procedure  requires  finer  control  for  increased  layer-to-layer 
thickness  uniformity. 

Cavity  Q  is  also  sensitive  to  absorption.  The  calculated 
spectrum  shown  in  Fig.  3  was  obtained  using  the  measured 

thicknesses  and  indices  and  an  approximated  complex  di- 
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FIG.  3.  Angle-dependent  reflectance  spectra  of  the  microcavity  resonator  for 
TE  and  TM  modes.  A  calculated  spectrum  is  shown  for  normal  incidence 
(dotted  line).  The  arrows  indicate  cavity  resonant  frequency. 


electric  constant  e".  The  fit  was  found  using  £*=0.04,  which 
is  reasonable  under  the  assumption  that  the  absorption  coef¬ 
ficient  a  =  (XQ)~l ,  where  X  is  the  cavity  resonant  fre¬ 
quency,  so  that  ^=2nTt0nyl02=2nTi02(aX/27r) 
=  (^Tiop/^0^O.O2. 

The  resonance  behavior  is  dependent  on  incident  angle. 
For  TE  polarization,  the  resonance  deepens  with  increasing 
incident  angle;  for  TM  polarization,  it  is  progressively  re¬ 
duced.  This  behavior  is  understood  in  light  of  the  reflectance 
bandwidths  for  the  interference  filter;  the  increasing  band¬ 
width  for  TE  modes  leads  to  improved  confinement  of  the 
resonance  at  higher  angles,  and  vice  versa  for  the  TM  modes. 
Increasing  the  angle  of  incidence  also  shifts  the  resonant  fre¬ 
quency  to  shorter  wavelengths.  The  resonance  tuning  may  be 
useful  for  modulation  of  signals.  The  cavity  resonance  is  also 
tuned  by  the  defect  layer  thickness.  Normal  incidence  spectra 
for  microcavities  with  different  cavity  thicknesses  are  shown 
in  Fig.  4. 

On  each  device,  surface  defects  such  as  flakes  and  cracks 
could  be  seen  via  optical  microscopy.  Regions  of  delamina¬ 
tion  both  from  the  substrate  and  from  the  uppermost  layers 
were  also  seen  by  XTEM.  These  defects  are  thought  to  be 
due  to  the  buildup  of  internal  stresses  induced  by  the  cyclic 
anneals  since  the  ratio  of  thermal  expansion  coefficients  is 
as\o2f aTio2^  1  *58  at  700  °C.]2  Additionally,  the  high  heating 
rate  corresponds  to  film  shrinkage  on  the  order  of  40%,° 
which  can  induce  significant  film  strain.  Improvements  to  the 
structural  integrity  of  the  final  devices  include  the  following: 
a  reduction  of  thermal  cycling  via  alternative  annealing  tech¬ 
niques,  and  the  inclusion  of  adhesion  or  accommodation  lay¬ 
ers  at  the  substrate  and  within  the  multilayer.  Another  effect 
of  uneven  surface  topology  is  its  contribution  towards  in¬ 
creasing  the  film  thickness  variation  for  subsequent  layers. 


1000  1500  2000  2500 

Wavelength  (nm) 


FIG.  4.  Normal  incidence  spectra  for  microcavities  with  different  defect 
layer  thicknesses  tdcfect .  The  arrows  indicate  cavity  resonant  frequency;  for 
1400  nm,  /dcfect=210  nm;  for  X=  1500  nm,  rdefcct=  270  nm. 


Clearly,  a  highly  controlled  deposition  and  annealing  process 
is  necessary  to  reduce  surface  roughness.  Despite  the  pres¬ 
ence  of  surface  defects,  our  sol-gel  interference  filters  are 
shown  to  be  sufficient  for  applications  that  require  high  re¬ 
flectance  over  a  large  surface  area. 

In  summary,  we  have  fabricated  via  the  sol-gel  method 
an  omnidirectional  reflector  using  Si02  and  Ti02  thin  films. 
Using  the  same  materials,  we  have  demonstrated  a  microcav¬ 
ity  with  resonance  at  X—1500  nm.  We  show  that  the  reso¬ 
nance  is  modulated  by  the  incident  angle  of  light  and  the 
defect  layer  thickness.  The  tolerance  of  reflectance  to  detun¬ 
ing  enables  the  sol-gel  method  to  be  used  for  fabrication  of 
omnidirectional  reflectors.  The  technique  must  be  highly 
controlled  for  fabrication  of  high  Q  microcavities. 
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Abstract — The  emergence  of  a  dielectric  omnidirectional  mul¬ 
tilayer  structure  [1]-J4]  opens  new  opportunities  for  low  loss 
broad-band  guiding  of  light  in  air*  We  demonstrate  the  ef¬ 
fectiveness  of  such  an  approach  by  fabricating  a  broad-band, 
low-loss  hollow  waveguide  in  the  10-pm  region  and  measuring  its 
transmission  around  a  90°  bend.  The  generality  of  the  solution 
enables  the  application  of  the  method  to  many  wavelengths  of 
interest  important  in  telecommunication  applications  as  well  as 
for  guiding  high-power  lasers  in  medical  and  other  fields  of  use. 

Index  Terms —  Dielectric,  high-power  lasers,  hollow  wavegu- 
sides,  light  conduits,  low-loss  broad-band  transmission,  medical 
lasers,  multimode  waveguide,  omnidirectional  refleetors,  optical 
fibers,  optical  confinement,  single-mode  waveguide. 


I.  Introduction 

GUIDING  light  in  dielectric  fibers  has  had  a  tremendous 
impact  on  many  aspects  of  our  life — we  rely  on  fiber 
optics  for  communications  as  well  as  for  illumination  and 
a  host  of  medical  applications.  The  typical  optical  fiber  has 
a  high  index  core  and  a  low  index  cladding  such  that  the 
light  is  confined  to  the  core  by  total  internal  reflection. 
Two  inherent  drawbacks  exist  in  this  approach:  the  first  is 
absorption.  Since  the  light  is  traveling  through  a  dense  medium 
for  long  distances,  material  absorption  becomes  significant 
even  in  low  loss  materials.  To  compensate  for  losses  the  fiber 
is  doped  with  erbium  which  is  used  to  amplify  the  signal. 
This  in  turn  limits  the  bandwidth  of  the  fiber  to  that  of  the 
narrow  erbium  excitation  lines.  The  other  weakness  follows 
from  the  confinement  mechanism — total  internal  reflection 
which  confines  light  only  of  a  limited  angle.  Conventional 
optical  fibers  cannot  guide  light  around  sharp  turns,  which  is 
especially  important  in  optical  integrated  circuits.  Light  guided 
in  a  hollow  waveguide  lined  with  an  omnidirectional  reflecting 
film  propagates  primarily  through  air  and  will  therefore  have 
substantially  lower  absorption  losses.  In  addition,  the  confine- 
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ment  mechanism  does  not  have  angular  dependence  allowing 
for  guiding  light  around  sharp  bends  with  little  or  no  leakage. 

Most  hollow  waveguides  fabricated  to  date  [5]-[7],  have 
internal  metallic  and  dielectric  layers.  It  has  been  shown  [8] 
that  the  addition  of  dielectric  layers  to  a  metallic  waveguide 
could  lower  the  losses  significantly.  In  contrast,  our  system 
is  an  all  dielectric  waveguide  which  confines  all  frequencies 
contained  in  its  omnidirectional  range.  In  principle  this  type  of 
structure  can  have  lower  losses  than  the  combined  metal  and 
dielectric  structure  since  the  waves  do  not  interact  with  a  lossy 
metallic  layer.  Although  our  proof  of  concept  demonstration 
involves  a  large  diameter  multimode  waveguide,  one  can 
fabricate  a  much  smaller  tube  that  could  in  principle  be  made 
to  support  a  single  mode. 

II.  Principle  of  Operation 

A  schematic  of  the  hollow  tube  is  presented  in  Fig.  I, 
as  well  as  the  index  of  refraction  profile.  In  a  realistic 
light  guiding  scenario  involving  many  bends  there  exist  no 
global  symmetries  and  thus  one  cannot  distinguish  between 
independent  TE  and  TM  modes.  Locally  one  can  define  a  plane 
of  incidence  with  respect  to  the  normal  to  the  film  surface  and 
the  incident  wave  vector.  Light  entering  into  such  a  tube  will 
invariably  hit  the  walls  many  times  and  explore  a  wide  range 
of  angle  of  incidence  of  both  polarizations  with  respect  to  any 
local  plane  of  incidence.  Since  the  air  region  is  bounded  by 
a  structure  that  has  a  gap  which  encompasses  all  angles  and 
polarizations  the  wave  will  be  reflected  back  into  the  tube  and 
will  propagate  along  the  hollow  core  as  long  as  kz  ^  0. 

III.  Sample  Preparation  Procedure 

A  Drummond  L92  mm  o.d.  silica  glass  capillary  tube  was 
cleaned  in  concentrated  sulfuric  acid.  The  first  tellurium  layer 
was  thermally  evaporated  using  a  LADD  30000  evaporator 
fitted  with  a  Sycon  Instruments  STM  100  film  thickness  moni¬ 
tor.  The  capillary  tube  was  axially  rotated  to  ensure  uniformity 
during  coating.  The  first  polymer  layer  was  deposited  by  dip 
coating  the  capillary  tube  in  a  solution  of  5.7  g  polystyrene 
DOW  61 5 APR  in  90  g  toluene.  The  next  layer  is  tellurium 
deposited  in  the  same  method  outlined  above.  The  subsequent 
polymer  layers  are  made  of  polyurethane  diluted  in  mineral 
spirits.  The  device  has  a  total  of  nine  layers,  five  Te  and  four 
polymer  and  a  total  length  of  10  cm.  The  layer  thickness 
are  approximately  0.8  /xm  for  the  tellurium  layer  (refractive 
index  4.6)  and  1.6  /im  for  the  polystyrene  layer  (refractive 
index  1.59).  An  optimal  design  will  vary  the  layer  thickness 
according  to  the  zeros  of  the  Bessel  functions.  Performance  as 
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Fig.  1.  Cross  section  of  the  hollow  waveguide  showing  the  hollow  core  and 
the  dielectric  films,  also  shown  is  the  index  of  refraction  profile  in  the  radial 
direction. 
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Fig.  2.  (a)  Measured  (dashed)  and  calculated  (dots)  normal  incidence  re¬ 
flectance  for  hollow  waveguide  in  the  radial  direction.  <b)  Calculated  grazing 
incidence  reflectance  for  the  TM  mode. 

well  as  the  layer  thickness  were  monitored  by  IR  spectroscopy. 
The  reflectivity  of  the  deposited  structure  was  measured  in 
the  radial  direction  using  a  Nicolet  FTIR  microscope  and  a 
variable  size  aperture,  to  ensure  domination  by  radial  reflec¬ 
tion.  The  coated  capillary  tube  was  then  inserted  in  a  heat 
shrink  tube  which  was  filled  with  silicone  rubber.  Finally,  the 
glass  tube  was  dissolved  using  concentrated  hydrofluoric  acid 
(48%).  The  resulting  hollow  tube  assembly  is  thus  lined  with 
the  mirror  coating  and  is  both  flexible  and  mechanically  stable. 


Fig.  3.  Hollow  tube  transmission  measurement  setup  on  the  spectrophotome¬ 
ter  (FTIR). 


Wavelength  (microns) 

Fig.  4.  Transmission  through  the  hollow  waveguide  around  a  90°  bend  as 
a  function  of  wavelength. 


layers.  In  addition  there  are  absorption  (8  /im)  peaks  due  to  the 
polyurethane.  Fig.  2(b)  is  the  calculation  of  the  reflectance  at 
grazing  incidence  for  the  TM  mode.  Since  the  omnidirectional 
frequency  range  is  defined  from  above  (high  frequency  edge) 
by  the  normal  incidence  gap  edge  (arrow)  and  from  below 
by  the  grazing  incidence  gap  edge  (arrow)  the  extent  of  the 
gap  is  completely  defined  by  these  two  data  points.  The  extent 
of  the  omnidirectional  range  for  the  parameters  used  in  this 
experiment  is  approximately  40%  [1],  [2]. 

The  transmission  through  the  tube  was  measured  using  a 
Nicolet  Magna  860  FTIR  bench  with  an  MCT/A  detector.  The 
transmission  was  measured  around  a  90°  bend  at  a  radius  of 
curvature  of  approximately  1  cm,  which  was  compared,  to 
the  straight  tube  transmission  to  correct  for  entrance  and  exit 
effects.  A  schematic  of  the  measurement  layout  is  presented 
in  Fig.  3. 

The  results  shown  in  Fig.  4  indicate  a  high  transmission 
around  the  90°  bend  for  a  spectral  band  that  corresponds  to 
the  omnidirectional  gap.  The  relatively  high  noise  level  in  the 
measurement  is  due  to  the  lack  of  purge.  This  measurement 
provides  a  proof-of-concept  indicating  the  low  loss  characteris¬ 
tics  and  guiding  abilities  of  the  all  dielectric  hollow  waveguide. 


IV.  Results  and  Discussion 

The  reflectance  measurements  and  simulations  are  shown 
for  normal  incidence  in  Fig.  2(a).  The  measured  gap  width  is 
smaller  than  predicted,  probably  due  to  microdefects  in  the  Te 
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An  All-Dielectric  Coaxial 
Waveguide 

M.  Ibanescu,1  Y.  Fink,2  S.  Fan,1  E.  L.  Thomas,2 
J.  D.  Joannopoulos1 

An  all-dielectric  coaxial  waveguide  that  can  overcome  problems  of  polarization 
rotation  and  pulse  broadening  in  the  transmission  of  optical  light  is  presented 
here.  It  consists  of  a  coaxial  waveguiding  region  with  a  low  index  of  refraction, 
bounded  by  two  cylindrical,  dielectric,  multilayer,  omnidirectional  reflecting 
mirrors.  The  waveguide  can  be  designed  to  support  a  single  mode  whose 
properties  are  very  similar  to  the  unique  transverse  electromagnetic  mode  of 
a  traditional  metallic  coaxiat  cable.  The  new  mode  has  radial  symmetry  and  a 
point  of  zero  dispersion.  Moreover,  because  the  light  is  not  confined  by  total 
internal  reflection,  the  waveguide  can  guide  light  around  very  sharp  corners. 


ilirough  the  empty  core  of  a  hollow  waveguide 
(7),  thus  minimizing  effects  associated  with 
material  nonlinearities  and  absorption  losses. 
Moreover,  because  confinement  is  provided  by 
the  presence  of  at  least  a  partial  photonic 
handgap,  this  ensures  that  light  should  be  able 
to  be  transmitted  around  a  bend  with  a  smaller 
radius  of  curvature  than  is  possible  with  the 
optical  fiber. 

Here  we  introduce  a  waveguide,  the  coax¬ 
ial  omniguide,  that  combines  some  of  the  best 
features  of  the  metallic  coaxial  cable  and  the 
dielectric  waveguides.  It  is  an  all-dielectric 
coaxial  waveguide  and  supports  a  fundamen¬ 
tal  mode  that  is  very  similar  to  the  TEM 
mode  of  the  metallic  coaxial  cable.  It  has 
a  radially  symmetric  electric  field  distribu¬ 
tion  so  that  the  polarization  is  maintained 


throughout  propagation.  It  can  be  designed  to 

Waveguides  are  the  backbone  of  modern  op-  optical  fiber  has  proven  to  be  undeniably  be  single-mode  over  a  wide  range  of  frequen- 

toelectronics  and  telecommunications  sys-  successful  in  many  ways,  it  is  nevertheless  cies.  In  addition,  the  mode  has  a  point  of 

terns.  There  are  currently  two  major,  and  very  plagued  by  two  fundamental  problems.  First,  intrinsic  zero  dispersion  around  which  a  pulse 

distinct,  types  of  waveguides  (metallic  and  because  the  fundamental  mode  in  the  fiber  can  retain  its  shape  during  propagation,  and 

dielectric)  that  are  used  in  two  separate  re-  has  an  electric  field  with  twofold  p-like  sym-  this  point  of  zero  dispersion  can  be  placed  in 

gimes  of  the  electromagnetic  spectrum.  For  metry,  the  polarization  of  light  coming  in  one  the  single-mode  frequency  window.  Finally, 

radio  frequencies,  the  metallic  coaxial  cable  end  of  the  fiber  and  the  polarization  coming  the  coaxial  omniguide  can  be  used  to  guide 

is  of  greatest  prominence  (/),  In  this  type  of  out  the  other  are  generally  completely  differ-  light  around  sharp  bends  whose  radius  of 

cable,  the  entire  electromagnetic  field  is  con-  ent.  This  leads  to  substantial  problems  when  curvature  can  be  as  small  as  the  wavelength 

fined  between  two  coaxial  metal  cylinders.  coupling  into  polarization-dependent  devices,  of  the  light.  This  waveguide  design  is  com- 

The  important  fundamental  electromagnetic  Second,  because  the  guiding  involves  total  pletely  general  and  holds  over  a  wide  range 

mode  of  a  coaxial  cable  is  the  transverse  internal  reflection,  it  is  not  possible  for  light  of  structural  parameters,  materials,  and  wave- 

electromagnetic  (TEM)  mode,  which  is  to  travel  in  a  fiber  with  a  sharp  bend  whose  lengths.  To  our  knowledge,  this  work  repre¬ 
unique  in  that  it  has  radial  symmetry  in  the  radius  of  curvature  is  less  than  3  mm  without  sents  the  first  successful  attempt  to  bridge  the 

electric  field  distribution  and  a  linear  rela-  significant  scattering  losses  (J).  For  light  at  disparate  modal  regimes  of  the  metallic  eo- 

tionship  between  frequency  and  wave  vector.  optical  wavelengths,  this  is  a  comparatively  axial  cable  and  the  dielectric  waveguides. 

This  gives  the  TEM  mode  two  exceptional  enormous  radius,  thus  limiting  the  scale  of  A  cross  section  of  a  traditional  metallic 
properties.  First,  the  radial  symmetiy  implies  possible  miniaturization.  coaxial  waveguide  is  shown  in  Fig.  I  A.  Light 

that  one  need  not  worry  about  possible  rota-  Recently,  however,  all-dielectric  wave-  is  confined  in  the  radial  direction  in  the  re- 
tions  of  the  polarization  of  the  field  after  it  guides  have  been  introduced  that  confine  opti-  gion  between  the  two  metal  cylinders  and 

passes  through  the  waveguide.  Second,  the  cal  light  by  means  of  one-dimensional  (ID)  (4)  travels  in  the  axial  direction  ( perpendicular  to 

linear  relationship  ensures  that  a  pulse  of  and  2D  (5,  6)  photonic  bandgaps.  Although  the  plane  of  the  figure).  In  a  simple  ray 

different  frequencies  will  retain  its  shape  as  it  single-mode  propagation  is  still  twofold  p-like  model,  propagation  of  light  through  the  co¬ 
propagates  along  the  waveguide.  The  crucial  symmetric,  these  new  designs  have  the  poten-  axial  cable  can  be  viewed  as  a  result  of 

disadvantage  of  a  coaxial  metallic  waveguide  tial  advantage  that  light  propagates  mainly  successive  specular  reflections  off  the  metal 

Is  that  it  is  useless  at  optical  wavelengths 
because  of  heavy  absorption  losses  in  the 
metal.  For  this  reason,  optical  waveguiding  is 
restricted  to  the  use  of  dielectric  materials. 

However,  because  of  the  differences  in 
boundary  conditions  of  the  electromagnetic 
fields  at  metal  and  dielectric  surfaces,  it  has 
not  previously  been  possible  to  recreate  a 
TEM-like  mode  with  all-dielectric  materials. 

Consequently,  optical  waveguiding  is 
done  with  the  traditional  index-guiding  (that 
is,  total  internal  reflection)  mechanism,  as 
exemplified  by  silica  and  chalcogenide  opti¬ 
cal  fibers.  Such  dielectric  waveguides  can 
achieve  very  low  losses  (2).  Although  the 


- - -  Fig.  1.  Schematics  of  coaxial  waveguide  cross  sections.  (A)  Traditional  metallic  coaxial  cable.  (B  and 

’Center  for  Materials  Science  and  Engineering  and  *-)  ^wo  embodiments,  of  type  A  and  B,  respectively,  of  the  all-dielectric  coaxial  omniguide  as 

Department  of  Physics,  Massachusetts  institute  of  described  in  the  text.  The  blue  cylindrical  dielectric  layers  correspond  to  high  dielectric  constant 

Technology,  Cambridge,  MA  02139,  USA.  ^Depart-  material,  and  the  green  cylindrical  dielectric  layers  correspond  to  low  dielectric  constant  material, 

ment  of  Materials  Science  and  Engineering,  Massa-  For  all  cases  (A  through  C),  light  can  be  confined  within  the  coaxial  region  (for  example,  r,  <  r  < 

chusetts  institute  of  Technology,  Cambridge,  MA  r0)  and  guided  along  the  axial  or  z  direction.  In  (A)  and  (B),  r \  =  3.00  a  and  rn  =  4.67  a,  whereas 

02139,  USA.  in  (C),  r%  =  0.40  a  and  r0  -  1.40  a. 
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walls.  The  dispersion  relations  for  the  first 
few  modes  supported  by  a  metallic  coaxial 
waveguide  are  shown  in  Fig.  2A.  For  defi¬ 
niteness,  the  inner  and  outer  radii  of  the 
waveguiding  region  are  taken  to  be  r.  =  3.00 
a  and  r0  =  4.67  a ,  respectively,  where  a  is  an 
arbitrary  unit  of  length  to  be  defined  later.  For 
any  value  of  the  wave  vector,  the  lowest 
frequency  mode  is  the  TEM  mode  for  which 
both  the  electric  and  magnetic  fields  are 
transverse  to  the  direction  of  propagation. 
This  mode  has  zero  angular  momentum, 
which  means  that  the  mode  is  invariant  under 
rotations  around  the  axial  direction.  Another 
useful  property  of  this  mode  is  its  constant 
group  velocity,  which  makes  it  dispersionless 
at  any  frequency.  The  other  modes  shown  in 
the  plot  are  transverse  electric  (TEm/ )  modes 
for  /  =  1  and  varying  angular  momenta  m  ( 8 ). 
The  cutoff  frequency  of  any  of  these  modes  is 
of  the  form 


where  /  is  the  solution  to  a  transcendental 
equation  for  each  value  of  the  angular  mo¬ 
mentum  m  and  for  each  polarization  (TE  or 

TM)  (i)* 

Designing  an  all-dielectric  waveguide 
with  similar  principles  of  operation  as  the 
metallic  coaxial  waveguide  is  not  straightfor- 


Reports 

ward,  because  the  boundary  conditions  at  a 
dielectric-dielectric  interface  differ  from 
those  at  an  air  (dielectric) -metal  interface.  In 
particular,  specular  reflections  cannot  be  ob¬ 
tained  on  a  dielectric-dielectric  interface 
when  the  ray  of  light  comes  from  the  region 
with  a  lower  index  of  refraction.  Thus,  it  has 
generally  been  assumed  that  an  all-dielectric 
coaxial  waveguide  cannot  be  designed  to  sup¬ 
port  a  TEM-like  mode,  even  in  principle. 
However,  recent  research  on  the  omnidirec¬ 
tional  dielectric  reflector  (9)  has  opened  new 
possibilities  for  reflecting,  confining,  and 
guiding  light  with  all-dieiectric  materials.  In¬ 
deed,  a  dielectric  hollow  waveguide  using 
this  principle  was  recently  fabricated  and 
tested  successfully  at  optical  wavelengths 
(4).  The  omnidirectional  dielectric  reflector, 
or  simply  the  dielectric  mirror,  is  a  periodic, 
multilayered  planar  structure  consisting  of 
alternating  layers  of  low  and  high  indices  of 
refraction.  This  structure  can  be  designed  so 
that  there  is  a  range  of  frequencies  at  which 
incoming  light  from  any  direction  and  of  any 
polarization  is  reflected.  Moreover,  the  elec¬ 
tric  fields  of  the  reflected  light  in  this  fre¬ 
quency  range  have  corresponding  phase 
shifts  that  are  quite  close  to  those  acquired 
upon  reflection  from  a  metal.  In  fact,  there  is 
a  frequency  for  each  angle  of  incidence  and 
each  polarization,  for  which  the  phase  shift  is 


identical  to  that  of  a  metal.  This  observation, 
together  with  the  fact  that  high  reflectivity  of 
the  omnidirectional  dielectric  mirror  is  main¬ 
tained  for  all  angles  of  incidence,  strongly 
suggests  exploration  of  the  possibilities  of 
using  an  omnidirectional  mirror  in  lieu  of  a 
metal  in  coaxial  cable  designs.  In  effect,  the 
omnidirectional  dielectric  mirror  provides  a 
new  mechanism  for  guiding  optical  and  in¬ 
frared  light  without  incurring  the  inherent 
losses  of  a  metal. 

In  a  cross  section  of  a  coaxial  waveguide 
(Fig.  IB),  the  metal  cylinders  from  the  left 
panel  have  been  replaced  with  cylindrical 
dielectric  layers  (10)  associated  with  an  om¬ 
nidirectional  mirror.  We  call  this  particular 
coaxial  waveguide  embodiment  coaxial  om¬ 
niguide  A.  The  parameters  of  the  dielectric 
layers  are  taken  from  the  hollow  waveguide 
experiment  (4)  to  be  the  following:  Layers 
shown  in  dark  blue  have  an  index  of  refrac¬ 
tion  /?,  =  4.6  and  a  thickness  d]  —  0.33  a, 
whereas  layers  shown  in  green  have  n2  =  1.6 
and  d2  —  0.67  a .  Here,  a  —  di  +  d2  is  the  unit 
length  of  periodicity  of  the  multilayered 
structure.  The  inner  (r.)  and  outer  (rj  radii  of 
the  coaxial  waveguiding  region  are  taken  to 
be  the  same  as  those  of  the  metallic  coaxial 
cable  described  earlier.  For  the  calculations 
presented  here,  we  have  set  the  index  of 
refraction  of  the  coaxial  waveguiding  region 


Traditional  Coaxial  Waveguide  Omnidirectional  Reflector  Coaxial  Omniquide-A 

n  so _ ... _ _ _  _ 


0  0.05  0.1  0.15  0.2  0.25  0.3  0  0.05  0.1  0.15  0.2  0.25  0.3  0  0.05  0.1  0.15  0.2  0.25  0.3 


Wavevector  kz  (2ir/a)  Wavevector  k2  (2it/a)  Wavevector  kz  (2n/a) 

Fig.  2.  Projected  band  structures  along  an  axial  direction.  (A)  Traditional  the  edge  of  the  light  cone.  The  horizontal  gray  fines  mark  the  boundaries 

metallic  coaxial  cable  with  inner  and  outer  coaxial  radii  of  r~  =  3.00  a  and  in  frequency  within  which  omnidirectional  reflectivity  is  possible.  (C) 

ro  =  4.67  a,  respectively.  The  red  bands  correspond  to  allowed  guided  Coaxial  omniguide  A  with  inner  and  outer  coaxial  radii  of  r}  =  3.00  a  and 

modes.  For  any  given  wavevector,  the  lowest  frequency  mode  is  a  TEM  ro  =  4.67  a,  respectively,  and  bliayers  consisting  of  indexes  of  refraction 

mode  characterized  by  a  perfectly  linear  dispersion  relation.  The  six  next  n1  =  4.6  and  nz  =  1.6  and  thickness  d1  =  0.33  a  and  d  =  0.67  a, 

highest  bands  correspond  to  transverse  electric  (TE,J  modes  with  /  =  1  respectively.  The  red  and  yellow  bands  Indicate  guided  modes  confined 

and  increasing  angular  momentum  m.  (8)  Omnidirectional,  reflecting,  to  the  coaxial  region  of  the  waveguide.  The  dashed  lines  indicate  modes 

all-dieiectric  multilayer  film.  Light-blue  regions  correspond  to  modes  for  with  less  than  20%  localization  within  the  coaxial  region.  There  is  dose 

which  light  is  allowed  to  propagate  within  the  dielectric  mirror,  and  correspondence  between  the  modes  labeled  m  =  1  to  m  =  6  and  those 

dark-blue  regions  correspond  to  modes  for  which  tight  is  forbidden  to  of  (A)  labeled  TEn  to  TE61.  Also,  the  yellow  m  =  0  mode  corresponds  to 

propagate  within  the  dielectric  mirror.  The  diagonal  black  line  identifies  a  TEM-like  mode,  as  discussed  in  the  text. 
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to  be  1 .  In  practice,  in  order  to  provide  struc¬ 
tural  support,  the  coaxial  waveguiding  region 
may  be  chosen  to  be  a  dielectric  with  a  low 
Index  of  refraction  without  affecting  the  fa¬ 
vorable  properties  of  the  coaxial  omniguide 

vt). 

Before  we  begin  our  investigation  of  the 
modes  supported  by  the  coaxial  omniguide, 
it  is  instructive  to  first  review  the  modes  of 
a  planar  omnidirectional  dielectric  mirror. 
The  projected  band  structure  of  the  omni¬ 
directional  mirror  is  shown  in  Fig.  2B.  The 
light  blue  regions  represent  allowed  propa¬ 
gation  modes  of  light  within  the  dielectric 
mirror.  The  dark  blue  region  represents 
modes  for  which  light  is  forbidden  to  prop¬ 
agate  within  the  dielectric  mirror.  The  thick 
black  line  identifies  the  edge  of  the  light 
cone,  and  the  horizontal  gray  lines  demar¬ 
cate  the  frequency  range  of  omnidirectional 
reflectivity.  It  is  precisely  within  this  range 
of  frequencies  that  one  would  expect  the 
coaxial  omniguide  to  support  modes  that 
are  most  reminiscent  of  those  of  the  metal¬ 
lic  coaxial  cable.  To  calculate  the  frequen¬ 
cies  and  field  patterns  of  the  modes  of 
coaxial  omniguide  A,  we  proceed  as  de¬ 
scribed  below. 

As  a  result  of  the  cylindrical  symmetry 
of  the  system,  there  are  two  good  “con¬ 
served  quantities”  that  can  be  used  to  spec¬ 
ify  and  classify  the  various  modes  support¬ 
ed  by  this  waveguide.  These  are  kxt  the 
axial  component  of  the  wave  vector,  and  m. 


the  angular  momentum  (jw  =  0,  1 ,  2  .  .  .). 
For  a  given  mode,  the  radial  and  angular 
components  of  the  electric  and  magnetic 
fields  can  be  calculated  from  the  corre¬ 
sponding  z  (axial)  components  (12).  For  a 
given  wave  vector  kz  and  angular  momen¬ 
tum  m,  the  axial  field  components  in  a  layer 
of  index  n  have  the  general  form 

,  F(z,r,<(>)  =  [A JJV)  +  BYm(£Tr)] 

X  (Cje***  +  **)&*-***  (I) 

where  F  stands  for  either  Ez  or  Hz;  Jm  and  Ym 
are  Bessel  functions  (13)  of  the  first  and  second 
kind,  respectively;  and  kT  is  a  transverse  wave 
vector  kT  =  V(n ta/c)2  —  k}. 

The  modes  of  the  coaxial  omniguide  are 
calculated  with  two  different  approaches.  The 
first  is  a  semianalytic  approach  based  on  the 
transfer  matrix  method  (14).  Starting  from 
Maxwell’s  equations,  the  z  components  of  the 
electric  and  magnetic  fields  in  each  layer  can 
be  written  in  the  general  form  given  by  Eq.  1. 
For  given  k2i  to,  and  m,  the  only  variables  that 
determine  the  E_  and  Hz  fields  are  the  four 
coefficients  in  front  of  the  Bessel  functions 
(two  for  Er  and  two  for  Hz),  The  boundary 
conditions  at  the  interfaces  between  adjacent 
layers  can  be  written  in  the  form  of  a  matrix 
equation 

lA\ 

B 
A' 


(2) 


Fig,  3.  Rower  density  in  the  electric  field  for  guided  modes  at  kx  =  0.19  (27r/a)  in  Fig.  2C,  (A  through 
D)  correspond  to  guided  modes  with  angular  momenta  m  =  0  to  m  -  3,  respectively.  The  color 
bar  indicates  that  power  increases  in  going  from  black  to  dark  red,  to  red,  to  orange,  to  yellow.  The 
blue  circles  identify  the  boundaries  between  the  various  dielectric  shells  and  are  included  as  a  guide 
to  the  eye.  Most  of  the  power  is  confined  to  the  coaxial  region  of  the  waveguide.  The  cylindrical 
symmetry  and  radial  dependence  of  the  m  —  0  mode  are  consistent  with  those  of  a  TEM  mode. 
[The  small  fluctuations  in  density  away  from  perfect  cylindrical  symmetry  in  (A)  are  a  consequence 
of  the  discrete  cubic  grid  used  in  the  computations.! 


where 


/  A  \ 
B 
A' 

\B'/, 


the  coefficients  that 


deter¬ 


mine  the  electric  and  magnetic  fields  in  the 
Jih  layer,  and  M  is  a  4-by-4  transfer  matrix 
that  depends  on  k_y  w,  m,  the  geometry  of  the 
layers,  and  their  Indices  of  refraction.  After 
calculating  the  electromagnetic  fields  for  a 
given  point  (&_,  u>),  we  find  the  resonant 
modes  by  examining  the  fractional  E-field 
power  confinement  in  the  coaxial  waveguid¬ 
ing  region;  that  Is 


dre(r)|E(r)|2w/ 


dre(r)|E(r)|2u> 


all  layers 


(3) 

The  second  approach  involves  a  numerical  so¬ 
lution  of  Maxwell’s  equations  in  the  frequency 
domain  with  the  use  of  the  conjugate  gradient 
method  within  the  supercell  approximation 
(15).  Supercells  of  size  (30  a  by  30  a  by  0.1  a) 
were  used,  leading  to  a  basis  set  of  about 
230,000  plane  waves.  Eigenvalues  were  consid¬ 
ered  converged  when  the  residue  was  less  than 
10~6.  The  results  of  both  approaches  were 
found  to  agree  to  better  than  0.1%. 

In  the  projected  band  structure  for  coaxial 
omniguide  A  (Fig.  2C),  the  red  and  yellow 
bands  represent  guided  modes  localized  with¬ 
in  the  region  defined  by  the  inner  and  outer 
coaxial  radii  of  the  waveguide.  The  dashed 
lines  represent  modes  with  less  than  20% 
localization  within  the  coaxial  region.  There 
is  close  correspondence  between  the  modes 
within  the  omnidirectional  reflectivity  range 
labeled  m  —  1  to  m  —  6  and  those  of  the 
coaxial  cable  labeled  TEn  to  TE61  (16).  The 
m  —  0  mode  appears  to  correspond  to  the 
TEM  mode.  Of  course,  for  a  coaxial  omni- 
gulde  with  a  limited  number  of  outer  shells, 
these  modes  can  only  exist  as  resonances. 
Nevertheless,  even  with  only  2.5  bilayers,  we 
find  that  they  can  be  extremely  well  localized 
resonances,  and  the  leakage  rate  decreases 
exponentially  with  the  number  of  shells.  The 
strong  localization  is  shown  in  Fig.  3,  Here, 
we  plot  the  power  density  in  the  electric  field 
for  the  four  lowest  frequency  modes  at  k.  = 
0. 19  (2 tt! a).  As  the  color  bar  indicates,  power 
Increases  in  going  from  black  to  dark  red,  to 
red,  to  orange,  to  yellow.  The  blue  circles 
identify  the  boundaries  between  the  various 
dielectric  shells  and  are  included  as  a  guide  to 
the  eye.  In  all  cases,  the  power  is  confined 
primarily  within  the  coaxial  region.  This  is 
particularly  true  for  the  m  ~  0  mode,  which  is 
also  cylindrically  symmetric,  just  like  the 
TEM  mode.  Although  it  is  well  known  that  a 
waveguide  consisting  only  of  dielectric  ma¬ 
terial  cannot  support  a  true  TEM  mode  (17% 
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we  find  that  the  m  =  0  mode  (which  is  a  pure 
TM  mode)  possesses  several  of  the  character¬ 
istics  of  the  TEM  mode.  First,  as  mentioned 
above,  it  has  zero  angular  momentum  and 
hence  a  radially  symmetric  electric  field  distri¬ 
bution.  Second,  the  electric  and  magnetic  fields 
within  the  coaxial  waveguiding  region  (where 
over  65%  of  the  power  is  concentrated)  are 
nearly  identical  to  those  of  the  metallic  coaxial 
cable;  lor  example,  the  predominant  compo¬ 
nents  are  Er  and  and  vary  as  l/r.  Finally,  at 
the  point  where  the  m  —  0  dispersion  curve 
(yellow  line)  crosses  the  light  line,  there  is  an 
exact  correspondence  between  the  electromag¬ 
netic  fields  of  the  coaxial  omniguide  and  the 
metallic  coaxial  cable,  inside  the  coaxial  region. 
Moreover,  the  derivative  of  the  group  velocity 
is  exactly  zero  near  this  point,  leading  to  nearly 
dispersionless  propagation  throughout  its  vicin¬ 
ity  (181 

The  characteristics  described  above  are  cer¬ 
tainly  the  attributes  one  would  hope  to  achieve 
in  order  to  overcome  problems  with  polariza¬ 
tion-rotation  and  pulse  broadening.  But  what 
about  single-mode  behavior?  The  bands  shown 
in  Fig.  2C  are  clearly  multimode;  that  is,  for  a 
given  frequency  there  are  two  or  more  guided 
modes  that  can  be  excited.  To  design  a  coaxial 
omniguide  that  can  support  single-mode  behav¬ 
ior,  we  need  only  readjust  our  structural  param¬ 
eters.  It  is  easier  to  keep  the  parameters  of  the 
bilayers  fixed,  so  that  the  omnidirectional  re- 


Wavevector  kz  (2rt/a) 

Fig.  4.  Projected  band  structure  for  coaxiat 
omniguide  B.  The  central  dielectric  rod  has  an 
index  of  refraction  n1  =  4.6  and  a  radius  r5  = 
0,40  a.  The  coaxial  waveguiding  region  has  an 
outer  radius  ro  =  1.40  a,  and  the  parameters  of 
the  outer  bifayers  are  the  same  as  those  used 
for  configuration  A.  The  yellow  dots  represent 
confinement  of  the  electric  field  power  of  more 
than  50%;  the  red  dots  represent  confinement 
between  20  and  50%;  and  the  dashed  tines 
represent  confinement  less  than  20%.  The  two 
white  boxes  identify  the  frequency  ranges 
where  the  m  =  0  band  exhibits  single-mode 
behavior. 


Reports 

flectivity  frequency  range  does  not  change. 
This  leaves  only  those  parameters  that  are  com¬ 
mon  to  both  the  coaxial  omniguide  and  the 
metallic  coaxial  cable:  the  inner  and  outer  radii 
of  the  coaxial  waveguiding  region.  Single¬ 
mode  operation  for  the  TEM-like  mode  will 
only  be  possible  if  all  other  modes  are  moved 
up  in  frequency  so  that  the  lowest  nonzero 
angular  mode  has  its  cutoff  frequency  inside  the 
bandgap.  To  do  this,  we  have  to  decrease  the 
inner  radius  of  the  coaxial  waveguiding  region. 
At  the  same  time,  the  thickness  of  the  bi layers, 
u,  should  remain  constant,  which  means  that  we 
can  no  longer  accommodate  three  bilayers  in 
the  inner  part  of  the  waveguide.  Actually,  the 
inner  radius  has  to  be  decreased  so  much  that 
we  are  forced  to  discard  the  periodic  structure 
in  the  inner  region  and  to  replace  it  with  a  single 
dielectric  rod.  Loss  of  the  inner-core  mirror 
structure  is  not  crucial,  however.  What  is  im¬ 
portant  is  to  add  a  thin  rod  of  dielectric  in  the 
core  in  order  to  avoid  the  l/r  divergence  of  the 
field  at  the  origin  and  to  use  a  dielectric  of  high 
enough  contrast  to  localize  the  TEM-like  mode 
in  the  coaxial  region.  This  approach,  however, 
will  not  work  if  r  >  a,  and  one  must  then  revert 
to  a  multilayer  core.  Testing  different  values  for 
the  inner  and  outer  radii  of  the  waveguiding 
region,  we  have  found  a  configuration  that  has 
the  desired  properties.  This  new  embodiment, 
coaxial  omniguide  B,  is  shown  in  Fig.  1C,  The 
central  dielectric  rod  has  an  index  of  refraction 
nx  =4,6  and  a  radius  rt  —  0.40  a.  The  coaxial 
waveguiding  region  has  an  outer  radius  rG  = 
1 ,40  a,  and  the  parameters  of  the  outer  bilayers 
are  the  same  as  those  used  for  configuration  A. 
In  configuration  B,  there  are  two  frequency 
ranges  where  the  waveguide  can  operate  in  a 
single-mode  fashion.  We  plot  the  dispersion 
curves  for  the  modes  supported  by  coaxial  Om¬ 
ni  guide  B  in  Fig.  4.  The  yellow  dots  indicate 
more  than  50%  confinement  of  the  electric  field 
power,  whereas  the  red  dots  represent  confine¬ 


ment  between  20  and  50%.  (The  dashed  lines 
indicate  confinement  that  is  less  than  20%.)  The 
two  white  boxes  identify  the  frequency  ranges 
where  the  m  =  0  band  is  single-mode.  A  com¬ 
parison  of  Figs.  2C  and  4  reveals  that  the  cutoff 
frequency  of  the  m  -  I  band  has  shifted  sig¬ 
nificantly  upward,  whereas  the  m  =  0  band 
remains  relatively  unchanged.  The  flatness  of 
the  m  =  1  band  {19)  enables  the  TEM-like  band 
to  be  single-mode  both  above  it  and  below  it  in 
frequency.  The  exact  values  of  the  parameters 
were  chosen  so  that,  in  the  middle  of  the  higher 
frequency  single-mode  window  [at  u>  =  0.205 
(2rrc/a)],  the  mode  is  also  dispersionless  (18). 

Figure  5  shows  the  distribution  of  the 
electric  field  components  for  the  m  =  0  mode 
of  coaxial  omniguide  B  at  k9  =  0.2  (2n7n) 
and  u>  =  0.203  (2n c/a).  Because  the  ( kz ,  ui) 
point  is  very  close  to  the  light  line,  the  elec¬ 
tric  field  in  the  waveguiding  region  is  almost 
completely  transverse  to  the  direction  of 
propagation.  (The  z  component  of  the  mag¬ 
netic  field  will  always  be  zero  because  this  is 
a  pure  TM  mode.)  The  field  distribution 
clearly  reveals  a  high  confinement  of  the 
mode  in  the  waveguiding  region,  as  desired. 
Moreover,  these  values  of  Ex  and  Ev  lead  to  a 
net  field  distribution  that  is  completely  radi¬ 
ally  symmetric,  consistent  with  an  angular 
component  that  is  exactly  zero.  All  the  fea¬ 
tures  mentioned  above  attest  to  the  close 
correspondence  between  the  m  =  0  mode  and 
a  pure  TEM  mode. 

There  are  several  additional  issues  associat¬ 
ed  with  the  coaxial  omniguide.  The  first  is  a 
practical  issue  involving  the  coupling  of  light 
into  the  coaxial  omniguide.  One  possible  meth¬ 
od  for  coupling  would  be  to  begin  with  an 
omniguide  with  a  very  thin  core  that  increases 
gradually  to  match  the  core  of  the  coaxial  con¬ 
figuration.  Because  the  electromagnetic  field  of 
a  laser  source  can  have  a  TEM^  mode,  this 
should  lead  to  efficient  coupling  into  the  TEM- 


Fig,  5.  Electric  field  for  the  mode  at  frequency  0,203  (2-ir c/a)  and  k2  -  0,2  (Z-nla)  in  Fig.  4.  (A 
through  C)  Electric  field  components  along  the  x,  y,  and  z  directions,  respectively.  The  color  bar 
indicates  that  large  positive  and  negative  values  are  shown  as  dark  red  and  dark  blue  regions, 
respectively,  whereas  white  areas  represent  regions  of  zero  values  of  the  electric  field,  and 
light-colored  areas  represent  regions  of  low  values  of  the  electric  field.  The  field  distribution  clearly 
reveals  a  high  confinement  of  the  mode  in  the  coaxial  waveguiding  region,  and  in  this  region  it  is 
nearly  completely  transverse  to  the  direction  of  propagation  (less  than  10“ 3  of  the  intensity  is 
along  z),  as  desired. 
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like  mode  of  the  coaxial  omniguide.  The  second 
issue  concerns  the  ability  to  transmit  informa¬ 
tion  across  a  broad  band  of  frequencies  with  the 
coaxial  omniguide.  This  is  determined  primar¬ 
ily  by  the  width  of  the  omnidirectional  reflec¬ 
tivity  range.  Another  issue  is  that  the  multitude 
of  adjustable  parameters  in  the  structure  of  a 
coaxial  omniguide  (the  index  of  refraction  and 
thickness  of  each  layer,  the  thickness  of  a  bi¬ 
layer,  the  waveguiding  region  inner  and  outer 
radii,  the  central  rod  index  of  refraction,  and  so 
on)  allows  for  great  flexibility  in  tuning  the 
waveguide  for  optimal  desired  performance 
(confinement  in  the  waveguiding  region,  width 
of  single-mode  window,  frequency  of  zero  dis¬ 
persion,  group  velocity,  and  so  on).  A  further 
important  issue  is  that  radial  confinement  of  the 
light  is  a  consequence  of  omnidirectional  re¬ 
flection  and  not  of  total  internal  reflection.  This 
means  that  the  coaxial  omniguide  can  be  used 
to  transmit  light  around  much  sharper  comers 
than  is  possible  with  the  optical  fiber.  Finally, 
the  radial  decay  of  the  electromagnetic  field  in 
the  coaxial  omniguide  is  much  faster  than  in  the 
case  of  the  optical  fiber  (with  only  1 0  bilayers, 
one  gets  a  decrease  of  electric  field  intensity  of 
about  six  orders  of  magnitude).  This  means 
that,  for  guided  light  of  the  same  wavelength, 
the  outer  diameter  of  the  coaxial  omniguide  can 
be  much  smaller  than  the  diameter  of  the  clad¬ 
ding  layer  of  the  optical  fiber  without  leading  to 
cross-talk  complications  in  waveguide  bundles. 
These  enabling  characteristics  of  a  substantially 
smaller  waveguide  bending  radius  and  smaller 
spacing  of  adjacent  waveguides  lead  to  the 
possibility  of  substantial  miniaturization  of  fu¬ 
ture  integrated  optical  devices  and  transmission 
lines. 
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Holes  in  a  Quantum  Spin  Liquid 

Guangyong  Xu,1  G.  Aeppli,2  M.  E.  Bisher,2  C.  Broholm,13* 

J,  F.  DiTusa,4  C.  D.  Frost,5  T.  Ito,6  K.  Oka,6  R.  L  Paul,3 
H.  Takagi,7  M.  M.  J.  Treaty2 

Magnetic  neutron  scattering  provides  evidence  for  nucleation  of  antiferromag¬ 
netic  droplets  around  impurities  in  a  doped  nickel  oxide- based  quantum  mag¬ 
net.  The  undoped  parent  compound  contains  a  spin  liquid  with  a  cooperative 
singlet  ground  state  and  a  gap  in  the  magnetic  excitation  spectrum.  Calcium 
doping  creates  excitations  below  the  gap  with  an  incommensurate  structure 
factor.  We  show  that  weakly  interacting  antiferromagnetic  droplets  with  a 
central  phase  shift  of  it  and  a  size  controlled  by  the  correlation  length  of  the 
quantum  liquid  can  account  for  the  data.  The  experiment  provides  a  quanti¬ 
tative  impression  of  the  magnetic  polarization  cloud  associated  with  holes  in 


a  doped  transition  metal  oxide. 

Spin  density  modulations  in  transition  metal 
oxides  are  receiving  huge  attention  because 
of  possible  connections  to  high-temperature 
superconductivity.  The  modulations  appear 
upon  introduction  of  charge  carriers,  through 
chemical  substitution,  into  an  insulating  and 
antiferromagnetic  parent  compound  and  tend 
to  be  static  when  the  carriers  are  frozen  and 
dynamic  when  they  are  mobile.  Evidence  for 
such  modulations  has  been  largely  confined 
to  materials  whose  magnetism  and  charge 
transport  are  quasi-two-dimensional  (1-4) 
and  whose  parent  insulators  are  ordered  anti- 
ferromagnets.  We  provide  evidence  for  anal¬ 
ogous  phenomena  in  a  quasi-one-dimenslon- 
al  oxide  (5),  Y2_^CaA.BaNi05,  for  which  the 
parent  is  a  spin  liquid  by  virtue  of  quantum 
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fluctuations  (6-10)  and  which  has  motivated 
considerable  theoretical  activity  (11-16), 

The  key  features  of  this  orthorhombic  ma¬ 
terial  are  the  chains  of  NiOs  octahedra  (Fig, 
l  A).  The  octahedra  are  comer-sharing,  which 
results  in  the  dominance  of  magnetism  (8) 
and  electrical  conduction  (5)  by  the  very 
simple  ..O-Ni-O-Ni-O.,  backbone.  The  mag¬ 
netic  degree  of  freedom  at  each  Ni  site  is 
the  spin  S  =  1  associated  with  the  3d8 
configuration  of  Ni2+,  Each  of  the  S  —  1 
ions  is  coupled  to  its  neighbors  through 
anti  ferromagnetic  (AFM)  superexchange 
through  shared  O2"  ions.  Replacing  the 
off-chain  Y3+  by  Ca2+  introduces  holes 
primarily  onto  apical  oxygen  atoms  in  the 
chains,  and  while  the  materials  remain  in¬ 
sulators,  doping  substantially  increases 
conductivity  at  finite  temperatures.  Thus, 
Y2_xCavBaNi05  is  a  one-dimensional  ana¬ 
log  of  the  cuprates,  where  off-Cu02 -plane 
chemical  impurities  donate  holes  to  the 
Cu02  planes. 

Magnetic  one-dimensionality  causes  the 
parent  compound  Y2BaNiOs  to  be  a  spin 
liquid  prevented  from  ordering  antiferromag- 
netically  by  quantum  fluctuations.  The  mate¬ 
rial  is  not  an  ordinary  paramagnet  with  heavi¬ 
ly  damped  spin  fluctuations,  but  rather  the 
magnetic  analog  of  superfiuid  4He  because  it 
has  a  macroscopically  coherent  quantum 
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